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INTRODUCTION 


The  Decision-Feedback  Equalizer  Simulation  (DFES)  program  is  written  in 
FORTRAN.  It  simulates  LMS  (Least  Means  Squares),  Kalman,  Rake,  and  fixed 
tap  versions  of  high  frequency  radio  channel  equalizers.  Considerable 
interest  in  this  program  has  been  shown  by  other  Government  laboratories, 
industry,  and  universities.  It  was  felt  necessary  to  preserve  the  following 
information  on  DFES  in  this  publication  in  order  to  ease  the  use  and/or 
modification  of  this  complex  program  by  future  users. 

The  material  herein  includes  basic  information  on  every  variable  in 
DFES  including: 

1.  Alphabetically  arranged  descriptions  of  all  variables. 

2.  Page  numbers  from  the  "Program  Performance  Specification" 
where  more  information  on  each  variable  may  be  found. 

3.  A  tabular  listing  of  all  variables  indicating  which  subprograms 
and  routines  set,  use,  or  output  each  variable. 

4.  A  listing  of  subroutines  and  functions  and  where  they  are  used 
in  the  program. 

This  document  is  useful  for  identifying  and  locating  variables  when  modifying 
and/or  debugging  the  program. 

In  the  following  sections,  DFES  refers  to  the  main  program  which  precedes 
the  initialize  parameters  subprogram. 


DFES  SUBPROGRAMS 


OFES 

Initialize  Parameters 
Update  Input 
Channel 
Interpolator 
Noise  Filter 
Forward  Filter 
Compressor 
Detector 

Differential  Decoder 
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FUNCTIONS 


ABS 

A I  MAG 

ALOG 

ALOG10 

AMOD 

CABS 

CEXP 

CMPLX 

CONJG 

COS 

DEXP 

ERFC 

EXIT 

EXP 

FLOAT 

IABS 

MOD 

RAN 

REAL 

SIGN 

SIN 

SINC 

SQRT 
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SUBROUTINES  -  WHERE  THEY  ARE  CALLED 


Subroutine 

Where  Cal  led 

BFILT 

detector 

DPGEN 

detector,  Sync,  Key 

FWATE 

detector 

GAURAN 

channel.  Noise 

KEY 

compressor,  Sync 

MAX 

Sync 

NOISE 

channel 

PIN1 

initialize  parameters 

SEMUL 

initialize  parameters 

SINC 

update  input 

SYNC 

noise  filter 

TAPER 

channel 
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Function 

ABS 

AIMAG 

ALOG 

ALOGI0 

AMOD 

CABS 

CEXP 

CMPLX 

CONJG 

COS 

DEXP 

ERFC 

EXIT 

EXP 

FLOAT 

IABS 

MOD 

RAN 

REAL 

SIGN 

SIN 

SINC 

SQRT 


FUNCTIONS  -  WHERE  THEY  ARE  USED 
Where  Used 

differential  decoder,  function  sine 

detector,  differential  decoder,  Fwate,  Bfilt,  Taper  (Rgen) 
Gauran 

detector,  differential  decoder,  Semul 

update  input,  detector,  differential  decoder,  Parin,  Fwate 

channel,  differential  decoder,  noise  filter,  forward  filter, 
detector,  Bfilt,  Fwate,  Sync,  Max 

channel,  noise  filter 

channel,  noise  filter,  detector.  Sync,  Fwate,  Bfilt,  Dpgen, 
Noise 

noise  filter,  detector,  differential  decoder,  Sync,  Fwate, 
Bfilt,  Key 

Gauran 

function  ERFC 
function  subroutine 

function  subroutine,  differential  decoder,  Sync,  Bfilt, 

Parin 

Semul 

forward  filter,  detector,  differential  decoder.  Taper 
(Read),  Fwate,  Semul,  Pinl,  Parin,  Sync 

Sync 

Sync,  Pinl,  Parin 
Gauran 

detector,  differential  decoder,  Fwate,  Bfilt 
detector,  Fwate,  Bfilt 
function  Sine,  Gauran 
function  Sine 

initialize  parameters,  channel,  detector.  Noise,  Sync, 
Gauran 
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DECISION  FEEDBACK  EQUALIZER  SIMULATION  (DFES)  PROGRAM  (FORTRAN) 


A 

local  variable  pages  (3-14,  3-45,  3-52,  3-67) 

Set  DFES 

set  &  used  detector 

output  differential  decoder  Set  &  used  Gauran 

Set  &  used  Sync  Set  &  used  Taper  (Rgen) 

A  =  Transmitted  PSK  digit  (Complex,  ARG' 

A  =  (1. ,  1. )  set  in  OFES 

The  transmitted  PSK  digit  is  then  formed  by  the  detector  as 
A  =  D 

A  =  D*EJ*A1  if  LTAPE  =  1 
Gauran  sets  A  to 

A  =  RAN  (NRAN1 ,  NRAN2) 

Taper  (Rgen)  sets 
A  =  D 


A 

array  pages  (3-99,  3-100) 

set  &  used  ERFC 
used  Max 

A  =  Complex  array  of  N  elements 


FHECEDING  PAfiS  BUNK-NOT  FliXSD 
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3. 


A1 


local  variable 
set  DFES 

used  noise  filter 
set  &  used  detector 
set  &  used  Sync 


pages  (4-1  back,  3-14,  3-65,  3-67, 
3-81) 

set  &  used  Sine 

set  &  used  Taper  (Rgen) 


Previous  PSK  Digit  at  Transmitter 

A1  is  the  complex  previous  transmitted  PSK  digit  after  encoding 
used  by  SNYC  to  differentially  encode  its  present  PSK  digit  A  obtained 
from  DPGEN  when  IDEC=1. 


A=A*A1*EJ 


where  EJ=(l+j)/2 

The  encoded  A  is  then  stored  in  A1  for  the  next  iteration.  When  the 
acquisition  decision  is  made,  the  latest  value  of  A1  is  stored  in  AHAT1 
as  the  previous  detected  PSK  digit  for  the  first  iteration  (NUM=1)  pass 
through  the  Detector  and  Differential  Decoder.  A1  is  updated  by  SYNC 
for  NA  bit  symbol  iterations. 

A1  is  set  in  DFES  to 

A1  =  (1.0,  1.0) 

Sine  sets  A1  to  The  detector  sets  A1  to  A1  =  A 

A1  =  ABS(PIX) 


4.  A2 

local  variable 
set  &  used  Sine 

A 2  is  set  in  Sine  as 

A2  =  ABS  (SIN  (PIX)) 


5.  AC 

local  variable 
set  &  used  channel 

AC  =  2(1-EC)EDC 
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6.  ADATA 


local  variable  pages  (3-82) 

set  &  used  Sync 

ADATA  is  the  modulation  PSK  symbol 
Sync  sets  ADATA  to 
ADATA  =  CONJG(A) 


7.  ADR 

local  variable 
set  &  used  Max 

ADR  is  set  in  MAX  as 

ADR  =  CABS  (A(MSET(I )))/CABS(A(MSET(l))) 


8.  AERR 

local  variable 
set  DFES 

used  differential  decoder 
AERR  =  .97723 


9.  AGCLB 

Common  block 
used  Pinl 

AGCLB  is  the  AGC  bandwidth  in  Hz,  real  default  =  10. 


10.  AGC LG 

local  variable 

set,  used,  &  output  Parin 

AGC LG  =0.01 
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11.  AHAT 


local  variable  pages  (4-1  back,  3-14,  3-46  to  3-50, 

set  DFES  3-52,  3-96,  3-97) 

set  &  used  detector 

used  &  output  differential  decoder 

used  Bfilt 

Detected  PSK  Digit 

AHAT  is  the  complex  detected  PSK  digit  developed  by  the  Detector 
each  bit  symbol  interation  after  the  acquisition  decision.  In  the 
Detector  it  is  used  to  compute  the  error  E.  The  DFES  main  program  sets 
AHAT=  1+j.  It  is  a  calling  sequence  argument  of  the  subroutine  BFILT 
where  it  is  used  to  form  ALPHA(I),  1=1,2,— NCB.  AHAT  has  the  values 
±1±  j . 


AHAT  is  set  by  DFES  to 
AHAT  =  (1.0,  1.0) 


12.  AHAT1 

Local  variable  pages  (3-14,  3-50) 

set  DFES 

set  noise  filter 

set  &  used  differential  decoder 

AHAT1  is  the  previous  AHAT  value 

The  noise  filter  sets  AHATl  to 

AHAT1  =  A1 

DFES  sets  AHATl  =  (1. ,  1.  ) 

The  differential  decoder  sets  AHATl  to 
AHATl  =  AHAT 


13.  AIX 

local  variable  pages  (3-38) 

set  &  used  noise  filter 

AIX  is  a  constant  for  the  2  pole  Butterworth  filter 
AIX  =  AN*CJ 

where  CJ  =  CEXP  (CMPLX  (0.,-CX)) 
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14.  AKC 

local  variable 
set  &  used  Fwate 

Fwate  sets  AKC  to 

AKC  =  1.  /  (1.  -  (6.28  *  KFLB  /  BSR)) 


15.  AL 

local  variable 
set  &  used  Semul 

AL  =  3. 


16.  ALGOR 

local  variable  pages  (3-55,  3-58,  3-92) 

set  &  used  initialize  parameters 

used  chanel  used  differential  decoder 

used  forward  filter  set  &  used  Parin 

used  detector  used  Fwate 

ALGOR  =  Algorithy  type  used  in  Forward  Filter  weight  adaptation 

(Integer  word  for  alphanumeric  input).  ALGOR  is  either  LMS 
(Least  Mean  Square),  KAL  (Kalman  algorithm),  or  FIX  (fixed 
weight  input).  Oefault  =  LMS. 

Initialize  parameters  sets  ALGOR  =  NALG 
Parin  sets  ALGOR  =  LMS,  default 


ALGOR 

NALG 

Weight  Adaptation  Method 

LMS 

1 

Least  Mean  Squares  algorithm 

KAL 

2 

Kalman  algorithm 

FIX 

3 

Weights  remain  fixed  to  initialization 
values 

RAKE 
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Rake  equalizer 


a)  local  variable  pages  (3-37) 

set  &  used  &  output  noise  filter 

7T  RDF 

ALPHA  =  —  -£75  =  normalized  filter  parameter 
V2  R5R 

b)  virtual  array  Previous  chip  values 

set  &  used  Sync  pages  (3-82,  3-96,  3-97) 

Sync  sets  ALPHA  as 

ALPHA  (I)  =  (0. ,0. ), 

ALPHA(l)  =  ALPHA  (I-NRB) 

ALPHA(I)  =  PN( I )*ADATA 

c)  array  ALPHA  (I)  =  AHAT*ALPHA  (I), 

1=1,2,. . .NCB. 

set  &  used  Bfilt 

ALPHA(I)  =  ALPHA  (>-NCB) 

Bfilt  sets  ALPHA  to  ALPHA(I)  =  CONJG(PN(I*NRC) 

ALPHA  (I)  =  (0.  ,0.  ) 

18.  AMAX 

local  variable 
set  &  used  Max 

AMAX  =  -1.0 

AMAX  =  CABS  (A(J)) 

19.  AN 

local  variable  pages  (3-37) 

set  &  used  noise  filter 

AN  =  2(1-ALPHA) 

20.  ARG 

local  variable 
set  &  used  Semul 

ARG  =  FLOAT  (I-l)/100. 

ARG  =  0/3  -  ALOG10(ARG) 


21.  ASTEP 


local  variable 
used  channel 
set  Pinl 

Pinl  sets  ASTEP  to 

ASTEP  =2.  *  PI  *  AGCLG/BSR 
and  if  NA>0  then  ASTEP  =  2.  *  PI/NA 


22.  B 

local  variable 
set  &  used  Gauran 

B  =  RAN  (NRAN1 ,  NRAN2) 


23.  BC 

local  variable 
set  &  used  channel 

BC  =  (1-2EC  +  2EC2 )  EDC2 


2A.  BDEL 

local  variable 
set  &  used  Pinl 

BDEL  as  set  in  Pinl 

BDEL  =  RSR  /  (3.  *  BRF)  +0.5 
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25.  BETA 


array  pages  (3-1  back,  3-17,  3-20, 

set  OFES  3-24,  3-29,  3-95  to  3-98) 

output  update  input  set  &  used  Bfilt 

output  differential  decoder  set  Parin 

OFES  sets  BETA(I)  =  0. 0,0.0) 

Backward  Filter  Tap  Array 

BETA  is  a  complex  array  dimensioned  BETA  (40).  For  NALG=1  or  2  it 
is  updated  each  bit  symbol  interation  following  the  acquisition  decision 
by  the  subroutine  BFILT  in  order  to  form  the  backward  filter  output  C. 
For  ALG0R=FIX,  the  BETA  remains  fixed  at  its  initial  values.  Within 
PARIN  the  array  is  indexed  as  BETA(I),  1=1,2. .. LTAP. 

BETA  is  a  calling  sequence  argument  of  both  PARIN  and  BFILT. 

Bfilt  sets  BETA  to 

BETA(I)  =  CMPLX  (RW.QW) 


26.  BFLB 

common  block  pages  (3-7  back,  3-55,  3-58 

used  Bfilt  3-95,  3-98) 

set,  used  &  output  Parin 

Backward  Filter  Loop  Bandwidth  in  Hz 

BFLB  is  an  input  parameter  to  PARIN  with  a  default  of  5.  It  is 
used  by  the  subroutine  BFILT  to  compute  BSTEP. 


27.  BIX 

local  variable  pages  (3-38) 

set  &  used  noise  filter 

BIX  is  a  constant  for  the  2  pole  Butterworth  filter 
BIX  =  BN  *  CJ 

where  CJ  =  CEXP  (CMPLX  (0. ,  -CX)) 
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28.  BKC 


local  variable 
set  &  used  Fwate 

Fwate  sets  BKC  to 

BKC  =  AKC  *  (l.-l./AKC) 


29.  BN 

local  variable 

set  &  used  noise  filter 

BN  =  2  ALPHA2 


30.  BRF 

common  block  pages  (5-1  back,  3-36,  3-37, 

used  noise  filter  3-53,  3-58) 

used  Pinl 

set,  used,  &  output  Parin 

RF  Bandwidth  in  Hz  of  2  Pole  Butterworth  Filter 


BRF  is  a  real  input  constant  to  PARIN  with  a  default  of  3840.  It 
is  used  by  the  Noise  Filter  subprogram  to  compute  the  constants  ALPHA 
and  FSPACE . 


31.  BRMS 

array  pages  (3-1  back,  3-13,  3-17,  3-21, 

set  OFES  3-26,  3-27,  3-54,  3-59) 

used  channel 

set,  used,  &  output  Parin 

Channel  Tap  Gain  Doppler  Spectrum  Standard  Deviation 

BRMS  is  a  real  array,  dimensioned  BRMS(4)  used  by  the  Channel 
subprogram  in  updating  the  Channel  Tap  Gain  array  H(I).  Indexing  is 
BRMS(I),  1=1,2 .CTAP.  Its  values  remain  fixed  to  the  initial  values 
input  by  the  subroutine  PARIN.  The  default  is  BRMS  (I)=0,  1=1 ,2, . . . CTAP. 
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common  block  pages  (3-5  back,  3-12,  3-15,  3-18, 

used  initialize  parameters  3-24,  3-25,  3-27,  3-29,  3-30, 

used  channel  3-38,  3-47,  3-55,  3-57,  3-58, 

used  noise  filter  3-66,  3-77,  3-87,  3-96,  3-98) 

used  differential  decoder 

used  Bfilt  used  Noise 

used  Fwate  used  Pinl 

Set,  used,  &  output  Parin 

Bit  Symbol  Rate  in  Hz 

BSR  is  a  real  input  constant  to  the  subroutine  Parin  with  a  default  of 
2400.  It  is  used  in  Parin  to  compute 

NTB  =  TSR/BSR 
NRB  =  RSR/BSR 
NCB  =  CSR/BSR 

Initialize  parameters  uses  it  to  set 

SIGMA  =  SQRT  (TSR/(BSR  *  2.  *  SNR)) 

The  channel  uses  BSR  to  compute 

DC  =  2.  *  PI  *  00P(J)/BSR 

EC  =  2.  *  PI  *  BRMS(J)/(BSR  *  SQ2) 

The  noise  filter  uses  it  to  set 

VSTEP  =  2  *  PI  *  NFLB  *  RSR  /  (BSR  *  BRF  *  NRB) 

The  differential  decoder  computes 

DRATE  =2.  *  BSR  *  (1.  -  FLOAT  (KADAPT)  /  KRST) 

Bfilt  uses  it  to  set 

BSTEP  =2.  *  PI  *  BFLB  /  BSR 

Fwate  uses  BSR  to  set 

AKC  =  1.  /  (1.  -  (6.28  *  KFLB  /  BSP)) 

Noise  uses  it  to  compute 

SVAR  =  TSR  /  (2.  *  BSR  *  SNR) 

Pinl  uses  BSR  to  compute 

DELTA  =  2.0  *  PI  *  FFLB  /  BSR 
ASTEP  =2.  *  PI  *  AGCLB  /  BSR 
ESTEP  =2.  *  PI  *  MSELB  /  BSR 
KDEL  =2.  *  PI  *  RGLB  /  BSR 


33.  BSTEP 


local  variable  pages  (3-97,  3-98) 

set,  used,  &  output  Bfilt 

BSTEP  =  (2n*BFLT)/BSR. 


34.  C 

local  variable  pages  (4-2  back,  3-45,  3-46,  3 

set  &  used  detector  3-96  to  3-98) 

output  differential  decoder 
set  &  used  Bfilt 

Backward  Filter  Output 

C  is  the  complex  backward  filter  output  argument  returned 
It  is  used  by  the  Detector  to  form  the- predecision  sample 

ZC=Z+C 

The  detector  sets  C  =  (0. ,0.  ) 

Bfilt  sets  C  as 
C  =  (0. ,0.  ) 

C  =  C  +  BETA(I)  *  CONJG  (GBACK(JSET(I))) 


35.  CAS 

Noise  filter  mode  indicator 
local  variable 

set  DFES  used  forward  filter 

used  noise  filter  set  &  output  Parin 

CAS  is  set  CAS  =  0  in  DFES  and  Parin 


36.  CC 

local  variable 
set  &  used  channel 


CC  =  2ECV2EC(1-EC)EDC 
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-52, 


by  BFILT 


37.  CIX 


local  variable  pages  (?-?8) 

set  &  used  noise  filter 

Constant  for  2  pole  Butterworth  filter 

CIX  =  CN  *  CJ  *  CJ 

where  CJ  =  CEXP  (CMPLX  (0. ,-CX) 


38.  CJ 

local  variable 

set  &  used  noise  filter 

CJ  is  set  in  the  noise  filter  as 

CJ  =  CEXP  (CMPLX  (0. ,-CX)), 

CJ  =  (0. ,0.), 

and  CJ  =  CJ  +  PN(K)  *  XFREQ 

39.  CLEAR 

local  variable 
set  &  used  Fwate 

CLEAR  =  0. 

CLEAR  =  1. 

if  CABS  (W( I ) )  >  0.2  then  CLEAR  =  0. 

40.  CMAG 

local  variable 

set  &  used  noise  filter 

CMAG  is  set  in  the  noise  filter  as 

CMAG  =  CABS  (V(I)) 

41.  CN 

local  variable  pages  (3-37) 

set  &  used  noise  filter 

CN  =  1-2  ALPHA  +  2  ALPHA2 
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42.  CON 

local  variable  pages  (3-89,  3-9G) 

set  &  used  Fwate 

NTAP  NLTAP  _ 

CON  =  (I  SBAR(i)  1  GBACK  (JSFT  (1-NTAP)))  *  VKAL(I)  *  BKC 
1=1  I=NTAP+1 

Fwate  sets  CON  to  CON  =  (0. ,0. ) 

if  I  <  NTAP  then 

CON  =  CON  +  CONJG  (SBAR(I))  *  VKAL(I)  *  BKC 
IF  I  >  NTAP  then 

CON  =  CON  +  CONJG  (GBACK  (JSET  (I-NTAP)))  *  VKAL  (I)  *  BKC 


43.  CSR 

common  block  pages  (3-6  back,  3-12,  3-55,  3-57, 

used  Taper  (Rgen)  3-58,  3-66) 

set,  used,  &  output  Parin 

Chip  Symbol  Rate  in  Hz 

CSR  is  a  real  input  constant  to  the  subroutine  PARIN.  It  has  a 
default  of  96K  and  must  be  an  integer  multiple  of  BSR. 


44.  CTAP 

common  block  pages  (3-17  back,  3-12,  3-15,  3-17, 

used  initialize  parameters  3-21,  3-23,  3-28,  3-54,  3-55, 

used  channel  3-58,  3-59,  3-60) 

set,  used,  &  output  Parin 

Number  of  Discrete  Channel  Paths 

CTAP  is  a  fixed  integer  input  to  the  PARIN  subroutine  with  a 
default  of  1.  It  has  the  range  1  <  CTAP  <  4.  If  the  input  value  of 
CTAP  is  greater  than  4,  PARIN  forces  CTAP=4.  The  channel  subprogram 
uses  CTAP  as  the  size  of  the  channel  arrays  H,  DELAY,  BRMS,  POW,  DOP, 
and  KSET.  The  Initialize  Parameters  subprogram  uses  it  to  set  the  size 
NQ  of  the  receiver  input  array  Q  as 

NQ=NTB+INT-NTR+KSET(CTAP) 
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45.  CVAR 

local  variable  pages  (4-2  back) 

set  channel 

Channel  Variance 


CVAR  is  the  variance  used  by  the  Channel  subprogram  as  an  argument 
to  the  random  number  generator  subroutine  GAURAN.  It  is  set  to  the  real 
va  1  ue 


CVAR=P0W( J )/2 


46.  CX 

local  variable 

set  &  used  noise  filter 

CX  is  set  as 

CX  =  2.  *  PI  *  1 1  *  F SPACE  /  RSR 
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47.  D 


h 


local  variable  pages  (4-2  back,  3-49,  3-50,  3-52, 

set  &  used  detector  3-67,  3-89) 

used  &  output  differential  decoder 

used  Taper  (Rgen) 

set  &  used  Dpgen 

Transmitted  PSK  Digit 

D  is  the  complex  transmitted  PSK  digit,  before  encoding,  returned 
to  the  Detector  by  DPGEN.  It  has  the  values  ±l±j.  The  Detector  uses  it 
to  set  the  transmitted  PSK  symbol  A=D.  The  Differential  Decoder  uses  it 
in  a  comparison  with  DHAT  to  update  the  error  counter  ERROR. 

The  following  process  is  performed  twice  to  generate  the  real  D(l) 
and  imaginary  D(2)  parts  of  the  complex  output. 

(a)  The  high  order  bit  I  is  extracted  from  MWORD. 

(b)  MWORD  is  shifted  one  place  to  the  left  by  extracting  the  low 
order  31  bits  and  multiplying  the  result  by  2. 

(c)  If  the  extracted  bit  I  is  a  1,  MWORD  is  replaced  by  an 
exclusive  OR  of  MWORD  and  JWORD. 

(d)  For  I  =  1,  D(J)  =  +1 
For  I  =  0,  D(J)  -  -1. 

The  detector  sets  D  to 
D  =  A  if  IDEC  =  0 

Dpgen  sets  D  to 
D(J)  =  MSIGN 


48.  DC 

local  variable 
set  &  used  channel 

DC  is  set  in  the  channel  as 

DC  =  2.  *  PI  *  DOP(J)  /  BSR 
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49.  DELAY 


array  pages  (3-21  back,  3-13,  3-17,  3-21, 

set,  used,  &  output  Parin  3-23,  3-28,  3-54,  3-59) 

Channel  Tap  Delays  in  Seconds 

DELAY  is  a  real  array  of  size  DELAY  (4)  used  by  the  channel  sub¬ 
program  in  updating  H(I),  1=1,2, — CTAP. 

Indexing  is  DELAY  (I),  I=1,2,--CTAP.  Its  values  remain  fixed  to 
those  input  by  the  subroutine  PARIN.  The  default  is 


DEUV  (I)  =  . 

I  =  1,  2,  --  CTAP 


50.  DELTA 

local  variable  pages  (5-1  back) 

output  update  input 
set  Pinl 
used  Fwate 

LMS  Algorithm  Step  Size 
Pinl  sets  DELTA  to 

DELTfl  =  2„|Et§  . 

It  is  the  calling  sequence  argument  for  the  LMS  algorithm  step  size  used 
by  FWATE  to  compute  the  weight  vector  W. 


51.  DERR 

local  variable 

set  &  used  differential  decoder 
DERR  =  ERROR 
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52.  DHAT 


Detected  source  digit 

local  variable  pages  (3-50) 

set  detector 

set,  used,  &  output  differential  decoder 
The  differentially  decoded  PSK  digit  is  given  by 
DHAT  =  AHAT  *  CONJG(EJ)  *  CONJG  (AHAT1) 

The  detector  defines  DHAT  as 
DHAT  =  AHAT 


53.  OOP 

array  pages  (3-21  back,  3-13,  3-17,  3-26, 

set  DFES  3-54,  3-60) 

used  channel 

set,  used,  &  output  Parin 
Doppler  Shift  in  Hz 

DOP  is  a  real  array  dimensioned  D0P(4)  used  by  the  channel  sub¬ 
program  in  updating  the  array  H.  Indexing  is  DOP(I),  I=1,2,--CTAP.  DOP 
is  a  fixed  real  parameter  input  to  PARIN,  with  a  default  of 

D0P( I )=0 ,  I=1,2,--CTAP. 


54.  DR 

local  variable  pages  (3-99) 

set  &  used  differential  decoder 
used  Max 

DR  is  the  Dynamic  Range  Threshold. 

The  differential  decoder  sets  DR  as 
DR  =  0 

if  ERROR  /  DERR  then  DR  =  DERR  -  ERROR 
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55.  ORATE 

local  variable  pages  (3-51,  3-52) 

set  &  output  differential  decoder 

ORATE  is  the  transmitted  data  rate  calculated  as 

DRATE  =  2  *  BSR  *  (1. -FLOAT(KADAPT)/KRST) 

by  the  differential  decoder. 


56.  E 

local  variable 
set  DEES 

set  &  used  detector 
used  Bfilt 

DEES  sets  E  to 

E  =  (0.0,  0.0). 

Adaption  Error  Sample 

E  is  the  complex  error  sample  computed  each  bit  symbol  iteration  by 
the  Detector  according  to  the  values  of  PSK  and  REF.  The  Detector  also 
uses  E  to  update  the  mean  square  error  MSE. 

The  detector  sets  E  to 

E  =  A-ZC 
E  =  AHAT-AC 
E  =  CMPLX(RE ,QE) 

E  =  E/2 

E  =  E  *  EMUl  ( IERR) 


57.  EBER 

local  variable  pages  (3-47,  3-48) 

set,  used  &  output  detector 

EBER  =  Estimated  bit  error  rate 

EBER  =  0. 5*ERFC(p) 


58.  EC 

local  variable 
set  &  used  channel 

EC  =  2n(BRMS(J))  /  V2BSR) 
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pages  (4-2  back,  3-13,  3-46  to  3-48, 
3-58,  3-95) 


59.  EOC 


local  variable 
set  &  used  channel 

EDC  =  CEXP  (CMPLX  (0. ,-OC)) 


60.  EDEL 

local  variable 
set  Sync 

EDEL  is  set  by  Sync  to 

EDEL  =  FLOAT  (INEXT  -  IMAX)  /  RSR 


61.  EJ 

local  variable  pages  (3-14,  3-50,  3-51,  3-66,  3-67) 

set  initialize  parameters 

used  detector 

used  differential  decoder 

set  &  used  Sync 

set  &  used  Key 

set  &  used  Taper  (Rgen) 

EJ  =  (1. ,1. )/2  in  initialize  parameters,  Sync,  Key,  and  Taper  (Rgen) 
EJ  =  (l+j)/2 


62.  EMSE 

local  variable 
set  &  used  Fwate 

EMSE  is  set  in  Fwate  to 

EMSE  =  1.  or  EMSE  -  0.1 


63.  EMUL 

array 
set  DFES 

set  &  output  Semul 

EMUL(l)  =1.0  Set  in  DFES 
In  Semul,  EMUL  is  set  to 
EMUL(l)  =  1. 

EMUL(I)  =  l.-SUM  *  EXP(-AL  *  ARG) 
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64.  ERFC 


local  variable 
set  Erfc 

Erfc  sets  ERFC  =  SUM 
The  function  computes 


ERFC(X)  =  —  J  e'y  dy,  X  >  0 
V n  X 

Using  the  series  approximation 


•  ^  I  -x^ 

ERFC(X)  =  1  A(I)*-r*e 

1=1 


65.  ERROR 


local  variable  pages  (3-1; 

set  DFES 

set,  used  &  output  differential  decoder 
Error  is  the  total  number  of  errors. 

ERROR  =  0.  Bit  error  counter  set  in  DFES 
The  differential  decoder  sets  ERROR  to 

ERROR  =  ERROR  +1  in  certain  cases 


pages  (3-13,  3-50,  3-51) 


66.  ESNR 


local  variable 

set,  used  &  output  detector 

Estimated  Signal  to  Noise  Ratio 


pages  (4-3  back,  3-47,  3-48) 


L_ 


ESNR  is  the  real  variable  for  the  estimated  SNR  computed  each  bit 
symbol  iteration  by  the  Detector  as 


ESNR  = 


2*MSE 


The  detector  sets  ESNR  to 


ESNR  =  -9.99E  +  32 


ESNR  =  20.  *  A LOG 10  (ESNR) 


ESTEP 

local  variable  pages  (3-15,  3-47) 

used  detector 
set  Pinl 

The  step  size  for  averaging  the  mean  square  error  is  computed 
ESTEP  =  2n(MSELB)/BSR. 


F 

local  variable  pages  (4-3  back,  3-36,  3-39,  3-47,  3-48) 

used  noise  filter 

set  &  used  detector 

output  differential  decoder 

Modified  Adaptation  Error  Sample 

The  Detector  computes  the  complex  error  F=E/NCB  to  be  used  as  a 
calling  sequence  argument  for  FWATE  and  BFILT  where  it  is  used  in  up¬ 
dating  W  and  BETA  respectively. 

The  detector  sets  F  to 

F  =  E 


FFLB 

common  block  pages  (3-8  back,  3-11,  3-55,  3-58) 

set,  used,  &  output  Parin 
used  Pinl 

Forward  Filter  Loop  Bandwidth  in  Hz 

FFLB  is  an  input  parameter  to  the  subroutine  PARIN  with  a  default 
value  of  5. 


FIX 

fixed  weight  input 

local  variable  pages  (3-3  back,  3-55) 

used  Parin 

algorithm 

no  numerical  value,  determines  algorithm 


T, 


71.  FSPACE 


local  variable  pages  (3-37) 

set,  used,  &  output  noise  filter 

FSPACE  =  BRF/2  =  Spacing  between  filter  center  frequencies 


72.  GBACK 

array  pages  (3-22  back,  3-88,  3-89,  3-96 

used  Fwate  to  3-98) 

set  Bfilt 


Backward  Filter  Signal  Arrav 


GBACK  is  a  complex  array  dimensioned  GBACK(40).  It  is  updated  each 
bit  symbol  iteration  with  the  latest  PSK  decision  and  chip  values  by  the 
subroutine  BFILT  which  uses  it  to  compute  the  backward  filter  output  C. 
The  subroutine  FWATE  uses  the  updated  GBACK  in  computing  the  variable 
CON  and  array  VKAL  in  the  Kalman  alogorithm  adaptation.  GBACK  is  a 
calling  sequence  argument  of  the  subroutines  BFILT  and  FWATE.  It  is 
indexed  as  GBACK(JSET(I)).  1=1,2,— LTAP. 


Bfilt  sets  GBACK  to 


GBACK  (JSET(I))  =  (0. ,0. ), 

GBACK(JSET(I))  =  GBACK  (JSET(I))  +  CONJG(PN( J*NRC) ) 
*  CONJG(ALPHA( J  +  JSET(I)-1))/NCB 


73.  GCON 

common  block 
set  DFES 

output  update  input 
set  &  used  channel 
used  interpolator 
output  detector 
set,  used,  &  output  Sync 
set  &  used  Parin 

Gain  Control  Constant 

When  the  bit  synchronization  subroutine  SYNC  is  required  for 
acquisition  (NA  >  0  and  SMODE  =  0),  GCON  is  computed  recursively  for 
IBS  =  NA  interations.  The  initial  value  of  GCON  is  1  set  in  DFES.  If 
NA  =  0,  which  forces  SMODE  =  1,  GCON  is  an  input  parameter  to  the  sub¬ 
routine  PARIN.  Whenever  SMODE  =  1,  the  vaiue  of  GCON  remains  fixed. 

The  channel  sets  GCON  to 

GCON  =  ( . . -ASTEP )*GCON+ASTEP/SQRT ( SPOW) 

Sync  sets 

GCON  =  GC0N/SQRT(SUM*2. ) 

Parin  sets 

GCON  =  I  if  GCON  =  0 


pages  (3-8  back,  3-11,  3-17,  3-19, 
3-29,  3-33,  3-56,  3-58) 


74.  H 

array  pages  (3-22  back,  3-12,  3-21,  3-26, 

set  DFES  3-27) 

set,  used,  &  output  channel 
set,  used,  &  output  Parin 

Channel  Tap  Gain 

H  is  a  complex  array  dimensioned  H(4)  for  the  channel  subprogram 
and  redefined  as  a  real  array  H(8)  in  the  parameter  input  subroutine 
PARIN.  Initial  values  are  input  by,  PARIN  with  the  default  H(l)=l,  all 
other  H(I)=0.  Indexing  is  H(I),  1-1,2, — CTAP.  The  array  H  is  updated 
each  bit  symbol  interation  for  ICHAN^O. 
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75.  HD 

Previous  channel  value 

array  pages  (3-?7) 

set  DFES 

set  &  used  channel 
set  Parin 

The  channel  subprogram  sets  HD  as 
HD(J)  =  YC 

and  the  initial  values  of  the  channel  are  stored,  i.e. 

HD( I )  =  H( I ) ,  1=1,2,. . -NCTAP2 

as  set  in  Parin  DFES  initially  sets  HD  to 

HD( I )  =  (0.0,  0,0) 
and  HD( 1)  =  (1.0,  0.0) 


76.  HP 

Preset  channel  value 

array 
set  DFES 

set  &  used  channel 
set  Parin 

Parin  sets  HP  to 

HP( I )  =  H( I )  where  I 

DFES  initially  sets  HP  to 

HP( I )  =  (0. 0,0.0) 

HP(1)  =  (1.0,  0.0) 


=  1.NCTAP2 

The  channel  sets 
HP( J)  =  H(J) 
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77.  I 


local  variable  pages  (3-38) 

used  DFES 

set  &  used  initialize  parameters 

set  &  used  differential  decoder 

set  &  used  update  input 

set  &  used  interpolator 

set  &  used  noise  filter 

set,  used,  &  output  forward  filter 

set,  used,  &  output  compressor 

set  &  used  Taper  (Read) 

set  &  used  Noise 

set  &  used  Erfc 

set  &  used  Parin 

set  &  used  Fwate 

set  &  used  Sync 

set  &  used  Bf i It 

set  &  used  Max 

set,  used,  &  output  Semul 

set  &  used  Dpgen 

set  &  used  Pinl 

example 

I  =  1,2,. ..MTAP. 


78.  IBDEL 

local  variable 

used  &  output  update  input 

used  forward  filter 

used  Fwate 

set  Pinl 

IBDEL  is  set  in  Pinl  to 
IBDEL  =  BDEL 


79.  IBETA 

algorithm 

local  variable  pages  (3-56,  3-59) 

set  &  used  Parin 

IBETA  =  Initial  BEAT  (I)  value  flag  (Integer). 
IBETA  =0,  BETA  (I)  =  default  values. 
IBETA=1 ,  BETA  (I)  from  hand  input. 
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80.  IBLOCK 


local  variable 

set,  used,  &  output  Taper  (Read) 

Taper  (Read)  sets  IBLOCK  to 

IBLOCK  =  0  and  IBLOCK  =  IBLOCK  +  1 


81.  IBRMS 

local  variable 
set  &  used  channel 

Channel  sets  IBRMS  to 

IBRMS  =  0  and  IBRMS  =  1 


82.  ICH 

array  pages  (4-3  back,  3-16,  3-65) 

set  Pinl 

used  Taper  (Rgen) 

Tape  Simulator  Chip  Shift  Register 

ICH  is  the  integer  shift  register  for  the  signal  tape  simulator 
subroutine  TAPER  (file  name  RGEN).  It  is  set  to  its  starter  value  by 
the  Initialize  Parameters  subprogram 

ICH  =  ICH1P 

TAPER-RGEN  uses  it  in  the  call  to  DPGEN  to  generate  the  chip 
values . 

Pinl  sets  ICH  to 

ICH(I)  =  ICHIP(I) 

ICH( 1+16)  =  ICHIP  (1+16) 


83.  ICHAN 

local  variable  pages  (3-14,  3-16) 

set  DFES 
used  channel 
#  used  Pari n 

If  ICHAN  =  1  indicating  a  non  ideal  channel  situation. 

If  the  default  condition  for  the  CHANNEL  subprogram  (3.4.3)  is 
detected  ICHAN  =  0. 
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84.  ICHIP 


array 

set,  used  &  output  Pinl 
used  Key 

Receiver  Chip  Shift  Register 

ICHIP  is  the  integer  shift  register  for  the  receiver  chip  data 
generation  of  the  PN  sequence. 

Pinl  sets  ICHIP  to 

ICHIP(I)  =  M0D( IL ,2) 

ICHIP( 1+16)  =  M0D(IH,2) 


85.  ICHP 

local  variable 
used  Taper  (Rgen) 


86.  IDATA 
array 

used  &  output  Taper  (Read) 


87.  IDEC 

common  block  pages  (3-20  back,  3-14,  3*15,  3-50, 

used  detector  3-66,  3*67) 

used  differential  decoder 

used  Sync 

used  Taper  (Rgen) 

set  Parin 

Differential  Decoder  Indicator 

Parin  sets  IDEC  =  1  for  REF=0,  and  IDEC=0  (no  decoding  required) 
for  REF  >  1.  When  IDEC=1,  the  error  rate  is  recalculated  as 

RATE=RATE/2 

If  IDEC=1,  the  Differential  Decoder  must  decode  to  obtain  DHAT  and 
the  SYNC  subroutine  must  encode  the  bit  symbol  returned  by  DPGEN. 
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local  variable  pages  (3-16) 

set  DFES 

Tape  Simulator  Previous  PSK  Symbol 
DFES  sets  IDEL1  to 
IDEL1  =  (1.0, 1.0) 


IDOP 

local  variable 
set  &  used  channel 

IDOP  is  set  in  channel  as 

IDOP  =  0  and  IDOP  =  1 


I  EOF 

local  variable  pages  (3-63,  3-64) 

set  &  used  Taper  (Read) 

I  EOF  =  0  End  of  file  indicator  returned  by 
Taper  (Read)  sets  IEOF  =  0 


I  EOT 

local  variable 
set  DFES 
used  channel 
se*  Noise 

IEOT  is  set 

I  EOT  =  0  in  DFES  and  NOISE 


TREAD 


=  1  for  end  of  NOISE(Tape) 


92.  I  ERR 


local  variable 
set  DFES 
used  detector 

set  &  used  differential  decoder 
set,  used,  &  output  Taper  (Read) 

IERR  =  1  in  DFES 

IERR  =  0  in  Taper  (Read) 

IERR  =  XERR  +  1  in  differential  detector 
and  if  IERR  >  50,  IERR  =  50 


93.  I F IX 

local  variable 
used  Parin 
integer  'FIX' 

no  numerical  value,  determines  algorithm 


94.  IH 

local  variable 
set  &  used  Pinl 

IH  =  24329  in  Pinl 

and 

IH  =  IH/2 


95.  II 

LOCAL  VARIABLE 

SET  &  used  noise  filter 

IT  is  set  in  the  noise  filter  as 

IT  =  I-NV 


96.  IK 

local  variable 
set  &  used  Fwate 

IK  is  set  in  Fwate  as 

IK  =  I  SET ( I )  +  K  +  IBDEL 
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97.  IL 

local  variable 
set  &  used  Pinl 

IL  is  set  in  Pinl  to 

IL  =  1432 

and 

IL  =  I L/2 


98.  I MAX 

local  variable 

set,  used,  &  output  Synr 

IMAX  is  set  to  IK *X  =  1 

and 

IMAX  =  MSVNC(l)  in  Sync 


99.  IMIN 

local  variable 
set  &  used  Sync 

IMIN  is  set  in  Sync  to 

IMIN  =  9999 


100.  INC 

local  variable 
set  &  used  Sync 

Sync  sets  INC  to 

INC  =  0, 

INC  =  INC  +  1, 

and 

INC  =  LM 


101.  INCH 

local  variable  pages  (3-56,  3-59) 

set  &  used  Parin 

INCH  =  Input  Channel  Flag  (Integer). 

0  =  Channel  default  parameter  values  used. 
1-Channel  parameter  values  from  hand  input. 
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102.  INEXT 


103. 


104. 


local  variable  pages  (3-85) 

set,  used,  &  output  Sync 

INEXT  is  set  by  Sync  to 

INEXT  =  1 

and 

INEXT  =  MSYNC(2) 


INT 

common  block  pages  (3-3  back,  3-12,  3-14,  3-15 

used  initialize  parameters  3-29,  3-54,  3-57,  3-79) 

used  update  input 
set,  used,  &  output  Parin 

Number  of  Interpolator  Samples 

INT  is  an  input  parameter  to  the  subroutine  PARIN.  It  is  an  odd 
integer  of  range  1  <  INT  <  11  with  a  default  of  5.  The  Initialize 
Parameter  subprogram  uses  it  to  set  the  Interpolator  half  span  interval 
INTI,  the  length  NQ  of  the  input  array  Q  and  the  length  NR  of  the 
receiver  sample  array  R  as  follows: 

INTI  =  (INT  -  l)/2 

NQ  =  NTB  +  INT  -  NTR  +  KSET  (CTAP) 

NR  =  NTB  +  INT  -  NTR 

In  the  bit  symbol  iteration  processing  it  is  used  to  set  the  number  of 
additional  input  samples  NRQ  to  be  provided  by  the  subroutine  TAPER  for 
IBS  =  1  to  NRQ  =  INT-NTR 

The  Initialize  Parameters  subprogram  will  reset  INT  for  INT  <  NTR 
to  NTR+1  for  NTR  even  or  NTR+2  for  NTR  odd.  From  that  point  on  the 
value  of  INT  remains  fixed. 

Parin  sets  INT  to 

INT  =  ((NTR+l)/2)  *2+1 


INTI 

common  block  pages  (3-13,  3-19,  3-33) 

set  initialize  parameters 
used  update  input 
used  interpolator 

INTI  =  (INT-l)/2  Interpolator  half  span 
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105.  INV 


local  variable  pages  (3-56,  3-59) 

set  &  used  Par  in 


INV  =  Input  V( I )  value  flag  (Integer). 
INV  =  0,  V( I )  =  default  values. 
INV  =  1,  V( I )  from  hand  input. 


*Li  ne  22:  INV 


0  =  default  V ( I )  values 
1  =  hand  input  of  V(I)  values 


This  line  is  input  only  for  MTAP  >  1 


106.  IP 

local  variable 

set  &  used  Taper  (Rgen) 

IP  is  set 

IP  =  ICHP  *  EJ  by  Taper(Rgen) 


107.  IPLUS 

local  variable 
set  &  used  Sync 

Sync  sets  IPLUS  =  0  and  IPLUS  =  1 


108.  IPOS 

local  variable 

set  &  used  Taper  (Read) 

IPOS  is  set  to  IPOS  =  0  and  IPOS  =  3  by  Taper(Read) 


109.  IPRIN 

local  variable  pages  (3-14,  3-22,  3-27,  3-40) 

set  update  input 

used  channel 

used  interpolator 

used  compressor 

used  noise  filter 

IPRIN  =  Print  parameter  from  UPDATE  INPUT  (3.4.2).  (Integer). 
IPRIN  =  0  Output  print  flag 
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110.  IQ 


111. 


112. 


local  variable 

set  &  used  interpolator 

The  interpolator  s6ts  IQ  as 

IQ  =  1  +  INTI  ♦  1 


IQSET 

local  variable 
set  initialize  parameters 
set  &  used  update  input 
set  noise  filter 

IQSET  is  set  IQSET  =  0  by  Update  Input. 

The  flag  IQSET  is  set  to  IQSET  =  1  in  the  initialize  parameters. 
IQSET=SMODE  as  set  in  the  noise  filter. 


I  RAKE 

local  variable 
used  Parin 


algorithm 


IRAKE  indicates  the  'RAKE'  a?gorithm~-RAKE  equalizer 
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113.  I  SET 


array 

used  initialize  parameters 
used  &  output  update  input 
used  &  output  forward  filter 
used  Fwate 

set,  used,  &  output  Sync 
set  &  used  Parin 

Transversal  Filter  Delay 

1SET  is  an  integer  array  of  non-negative  values  dimensioned 
ISET(IOO).  It  may  be  input  by  PARIN  or  computed  by  the  subroutine  SYNC, 
after  which  its  values  remain  fixed.  ISET(l)  defines  the  main  tap  of 
the  forward  filter  and  is  used  by  the  Initialize  Parameters  subprogram 
to  compute  the  parameter 

NSHIFT  =  NTR*ISET(1) 

Indexing  is  ISET(I),  1=1,2, - NTAP.  In  the  Forward  Filter  subprogram 

I  SET ( I )  is  used  in  the  index  value  of  X  for  computing  the  forward  filter 
output  Y.  The  subroutine  FWATE  also  uses  I SET ( I )  in  index  value  of  X 
for  updating  the  array  SBAR.  ISET  is  a  calling  sequence  argument  of  the 
subroutines  SYNC,  FWATE,  and  PARIN. 

ISET ( 1)  defines  main  tap  of  forward  filter. 


pages  (3-23  back,  3-13,  3-16,  3-18 

3-20,  3-41  to  3-43,  3-53,  3-81, 
3-83  to  3-85,  3-87,  3-89) 


114.  ISMAX 

local  variable  pages  (3-18,  3-42) 

set  &  used  update  input 

The  span  of  the  forward  filter  is  computed  as 

ISMAX  =  MjX  ISET(I) 

Update  input  sets  ISMAX  to 

ISMAX  =  0 

ISMAX  =  ISET(I)  if  ISET  (I)  >  ISMAX 


115.  I  SNR 

local  variable 
set  DFES 
used  Parin 

ISNR  is  set  in  DFES  to  the  values  SNR  is  to  be  set  at. 
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116.  IT 


local  variable 

set  initialize  parameters 

used  Semul 

Initialize  Parameters  sets  IT  as 
IT  =  KFLB 


117.  ITRY 

local  variable 

set  &  used  Taper  (Read) 

set  &  used  Parin 

ITRY  is  set  by  Taper  (Read)  as 

ITRY  =  0  and  ITRY  =  ITRY  +  1 

Parin  sets  ITRY  to 

ITRY  =  0  and  ITRY  =  1 


118.  IW 

local  variable 
set  &  used  Parin 

IW  =  Initial  Weight  Value  flag  (Integer) 

IW  =  0,  W( I )  =  default  values. 

IW  =  1,  W(I)  initial  values  from  hand  input 
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119.  IX 


local  variable  pages  (3-18,  3-36,  3-37) 

set  and  output  initialize  parameters 
used  noise  filter 


The  number  of  bit  symbol  iterations  to  be  executed  with  no  signal 
present  is  computed. 


If  NOSIG  =  0  (Signal  Present),  IX  : 
If  NOSIG  =  1  (Signal  Absent),  IX  = 


0 

3( BSR ) 
2n(NFLB) 


The  simulation  automatically  cycles  through 


IX  =  [3*RSR/(2n*NFLB*NRB)] 

x  <  [x]  <  x  +  1,  [x]  integer 

The  NOISE  Filter  also  has  the  number  of  adaptation  cycles  input  from 


UPDATE  INPUT  (3.4.2); 

IX  =  3*RSR/(2. *P1*NFLB*NRB) 


120.  IXF 

local  variable 
set  nose  filter 


IFX  =  J  +  1 
IFX  =  J 


set  &  used  Fwate 
set  &  used  Bf i It 
set  &  used  Max 
set  &  used  Semul 
set  &  used  Dpgen 
set  &  used  Parn 

set  &  used  Sync 

example:  J  =  1,  CTAP 


The  noise  filter  sets  IXF  as 
IXF  =  1 

IXF  =  (NRB  +  K-l)  *  NRB  +  I 
IXF  =  (K-l)  *  NRB  +  1 


121.  J 

local  variable 

set  &  used  update  input 

set  &  used  channel 

set  &  used  noise  filter 

set,  used,  &  output  forward  filter 

set  &  used  Taper  (Read) 
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122.  JBLOCK 


123. 


124. 


125. 


126. 


local  variable 

set,  used,  &  output  Taper  (Read) 

JBLOCK  is  set  by  Taper  (Read)  as 

JBLOCK  =  0  and  JBLOCK  =  JBLOCK  +  1 


JBS 

local  variable 
set  &  used  channel 

Channel  sets  JBS  to 

JBS  -  RJBS 


JFACT 

local  variable 
set  &  used  Semul 

JFACT  is  set  by  Semul  as 

JFACT  =  1  and  JFACT  =  JAFACT  *  (J) 


JH 

local  variable 
set  &  used  Pinl 

Pin  1  sets  JH  as 

JH  =  "7702  and  JH  =  JH/2 


JL 

local  variable 
set  &  used  Pinl 

Pin  1  sets  JL  as 

JL  =  "27607  and  JL  =  JL/2 
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127.  JMAX 


local  variable 
set  &  used  Bf i 1 t 

Bfilt  sets  JMAX  as 

JMAX  =  0 

JMAX  =  JSET(l) 

JMAX  =  JMAX  +  NCB-1 

JMAX  =  ( JMAX/NCB)  *  NCB  +  NCB 


128.  JP 

local  variable 
set  &  used  Semul 

JP  is  set  by  Semul  as 

JP  =  1,  M+l 


129.  JRAN1 

Integer  starter  for  Gauran  used  by  noise 

local  variable 
set  DFES 

DFES  sets  JRAN1  =  0 


130.  JRAN2 

Integer  starter  for  Gauran  used  by  noise 

local  variable 
set  DFES 

DFES  sets  JRAN2  =  0 
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131.  JSET 


array  pages  (3-23  back,  3-53,  3-59,  3-81, 

set  DFES  3-85,  3-88,  3-89,  3-95,  3-97) 

set,  used,  &  output  Sync 

used  Fwate 

used  &  output  Bfilt 

set  Parin 

Backward  Filter  Delay 

JSET  Is  an  integer  array  of  non-negative  values  dimensioned 
JSET(120).  It  may  be  input  by  PARIN  or  computed  by  SYNC.  Indexing  is 
JSET(I),  1=1,2, — LTAP .  JSET  is  a  calling  sequence  argument  to  the 
subroutines  SYNC,  BFILT,  FWATE,  and  PARIN.  FWATE  uses  JSET  as  an  index 
for  the  array  GBACK  and  BFILT  uses  it  in  the  index  for  the  GBACK  and 
ALPHA  arrays. 


132.  JTIME 

local  variable 
set  DFES 

used  differential  decoder 

used  Parin 

used  update  input 

JTIME  is  set  in  DFES  to  values  desired;  such  as,  JTIME  =  1,3 
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133.  JWORD 


a)  array  pages  (3-69) 

used  Dpgen 

set  &  output  Pinl 

Pinl  sts  JWORD  as 

JWORD(I)  M0D( JL , 2 ) 

JWORD( 1+16)  =  M0D(JH,2) 

b)  common  block  pages  (3-14  back) 

set  &  output  Pinl 

used  Dpgen 

Polynomial  for  the  Receiver  Message  and  Chip  Random  Number 
Generator 

JW0RD=I*(2**16)+J  Pinl  sets  JWORD  to 

where  JWORD(I)  =  M0D(JL,2) 

I  =  7702g  JW0RD( 1+16)  =  M0D(JH,2) 

J  =  27607g 

JWORD  has  the  actual  representation  of  1760427607 

The  message  and  chip  random  number  generator  DPGEN  uses  JWORD 
in  an  exclusive  OR  with  the  message  in  chip  shift  register  MWORD  to 
produce  a  new  shift  register  value  when  the  last  shift  resulted  in 
a  carry  of  "1. " 

The  integer  JWORD  is  the  COMMON  block  32  bit  polynomial  used 
for  both  message  and  chip  data  generations. 

JWORD  =  1760427607  (base  8) 


134.  K 

local  variable  pages  (3-38,  3-40,  3-41,  3-43,  3-64) 

set  &  used  channel 

set  &  used  interpolator 

set,  used,  &  output  noise  filter  set  &  used  Sync 

set,  used,  &  output  forward  filter  set  &  used  Fwate 

set  &  used  Taper  (Read)  set  &  used  Dpgen 

examples  K  =  NRB ,  NRB-1,  ...  2,1 

the  last  K  unused  data  samples  (K=HIGH-NVIN) 
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135.  KADAPT 


136. 


137. 


common  block  pages  (3-20  back,  3-11,  3-15, 

used  forward  filter  3-49,  3-56) 

used  detector 

used  differential  decoder 

used  Pinl 

set,  used,  &  output  Parin 

Number  of  Adaptation  Cycles  for  Kalman  Algorithm 

KADAPT  is  a  fixed  integer  parameter  input  to  PARIN  when  ALG0R=KAL. 
It  has  a  default  of  1000.  KADAPT  is  used  to  compute  RATE  and  DRATE. 

Parin  uses  it  to  set  KEND=KADAPT+KVAR. 


KAL 

al gori thm 

local  variable 
used  Parin 


pages  (3-55) 


KAL  (Kalman  algorithm) 

No  numerical  value,  determines  algorithm. 


KALA 

local  variable  pages  (3-16) 

set  DFES 

set  &  used  initialize  parameters 
used  channel 

DFES  sets  KALA  =  0 

Initialize  parameters  sets  KALA  =  2  *  NTAP 
if  NALG  =  2  and  SMODE  =  1 

This  is  accomplished  with  the  flag  KALA  which  is  used  to  fix  the 
weights  (NALG=3)  for  this  period  of  time. 
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138.  KEND 


common  block  pages  (3-19  back,  3-15,  3-91) 

used  Fwate 

set  &  used  Parin 

Kalman  Adaptation  Cycle  Limit 

Subroutine  PARIN  sets  the  integer  constant 

KEND-KADAPT+KVAR 

It  is  used  by  FWATE  in  a  decision  making  comparision  with  KNUM. 


139.  KF  LB 

Kalman  filter  loop  bandwidth.  Default  =  0 

common  block 
set  DFES 

used  initialize  parameters  used  Fwate 

used  update  input  set  &  output  Parin 

used  detector 

DFES  sets  KFLB  =0.0 

Parin  also  sets  KFLB  =  0 


140.  KKN 

common  block 

sets  DFES 

set  &  used  Fwate 

KKN  is  set  by  DFES  as 

KKN  =  1 

KKN  is  set  by  Fwate  as 
KKN  =  2 
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141.  KLMS 


common  block  pages  (3-19  back,  3-88,  3-96,  3-98) 

set  detector 
used  Bfilt 

set,  used,  &  output  Fwate 
Kalman/LMS  Adaptation  Indicator 

KLMS  is  an  integer  variable  with  a  value  of  0  or  1.  Initially  when 
NUM= 1 ,  FWATE  sets  KLMS=1  for  NALG=1  and  KLMS=0  for  NALG=2.  When  KLMS=1, 
FWATE  updates  the  weight  vector  W  using  the  least  mean  squares  algorithm 
For  KLMS=0,  the  Kalman  algorithm  is  used.  For  NALG=2  FWATE  leaves  KLMS^O 
until  KNUM  >  KEND  or  KNUM=KRST  when  it  sets  KLMS=1.  For  each 
KVAR  <  KNUM  <  KEND,  KLMS=0. 

In  the  subroutine  BFILT  the  method  used  to  compute  the  array  BETA 
depends  on  the  value  of  KLMS. 


142 .  KNUM 

common  block  pages  (3-18,  3-19  back,  3-14, 

set  DFES  3-90,  3-91) 

set,  used,  &  output  Fwate 

Kalman  Adaptation  Bit  Symbol  Iteration  Counter 

The  integer  variable  KNUM  is  set  to  0  by  the  DFES  main  program. 
When  NALG=2,  the  subroutine  FWATE  increments  KNUM  by  1  each  bit  symbol 
iteration  until  KNUM=KRST.  At  this  point,  it  is  reset 

KNUM=0 

FWATE  sets 

REF=3  for  KNUM=KVAR 

REF=0  for  KNUM=KEND 

KLMS=0  for  KNUM  <  KEND 

or  KNUM  >  KRST 

KLMS=0  for  KNUM  >  KVAR 

KNUM  is  used  computing  the  FWATE  variables 

AKC=F  LOAT ( KNUM+1 )/KNUM 

BKC=1. /KNUM 
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143.  KRS  T 


144. 


common  tlock  pages  (3-19,  3-20  back,  3-49,  3-56, 

used  detector  3-58,  3-88,  3-91) 

used  differential  encoder 
used  Pinl 

set,  used,  &  output  Parin 
used  Fwate 

Number  of  Cycles  for  Kalman  Algorithm  to  Restart 

KRST  is  a  fixed  integer  input  to  PARIN  when  ALG0R=KAL.  It  has  a 
default  value  of  100.  KRST  is  used  in  computing  the  error  rate  RATE  and 
actual  data  rate  DRATE  when  NALG=2.  FWATE  uses  KRST  in  a  comparison 
with  KNUM  in  order  to  set  KLMS. 


KSET 

array  pages  (3-11,  3-15,  3-17,  3-23,  3-24) 

set  DFES 

used  initialize  parameters 

used  channel 

set  &  output  Parin 

DFES  sets  KSET  to 

KSET ( I )-0 

If  ICHANjfO,  the  channel  transversal  filter  delays  are  computed  as 

KSET ( I )  =  TSR*DELAY( I )+0. 5  1-1,2, ... CTAP 

KSET(I)  =  Transversal  filter  delay  of  Ith  channel  tap  expressed  in  units 
of  1/TSR.  (Integer) 

KSET(I)>0=1,2,. . .CTAP. 

If  a  nontrivial  channel  is  selected  ICHAN=1  and  the  subprogram 
computes  the  set  of  integers  corresponding  to  the  number  of  tape 
sampling  intervals  for  each  path  delay,  i.e., 

KSET ( I )  =  TSR*DELAY( I )  +0.5.  =  XKSET 

as  in  Parin 
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145.  KSTEP 


local  variable  pagu-  (4-4  back,  3-89  to  3-92, 

set  &  used  Fwate  3-96,  3-98) 

used  Bfi It 

Kalman  Step  Variable 


KSTEP  is  a  real  variable  computed  by  FWATE  for  use  in  its  Kalman 
algorithm  adaptatior  of  the  w.  ight  vector  W  and  in  the  Kalman  algorithm 
update  of  the  BETA  a-'-ay  in  BFILT.  KSTEP  is  a  calling  sequence  argument 
of  FWATE  and  BFILT. 

Fwate  sets  KSTEP  to 

KSTEP  =  BKC/(EMSE  +  REAL(CON) ) 


146.  KSYNC 

local  variable  pages  (3-18) 

set  initialize  parameters 
used  update  input 
set  noise  filter 

Synchronization  flag  is 

KSYNC  =  1. 


147.  KVAR 

common  block  pages  (3-20  back,  3-11,  3-15,  3-56, 

set,  used,  &  output  Parin  3-58,  3-88,  3-90,  3-91) 

used  Fwate 

Number  of  Iterations  for  Adaptation  of  Kalman  Inverse 
KVAR  is  a  fixed  integer  parameter  input  to  PARIN. 

Parin  sets  KVAR  =  0 
PARIN  uses  it  to  set 
KEND=KADAPT +KVAR 

The  subroutine  FWATE  uses  it  to  set 
REF=3  for  KNUM=KVAR 

and 

KLMS=0  for  KNUM  >  KVAR 
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148.  LIMIT 


local  variable  pages  (3-82) 

set  &  used  Sync 

LIMIT  =  NPN  +  NRB-1 


149.  LM 

local  variable 
set  &  used  Sync 

Sync  sets  LM=0  and  in  special  cases  to  LM=LM+1 

150.  LMS 

algorithm 

local  variable  pages  (3-55) 

used  Parin 
no  numerical  value 
determines  algorithm 

LMS  (Least  Mean  Square) 

Default  =  LMS 

151.  LOUT 

common  block 
set  DFES 

LOUT  is  set  in  DFES  as  LOUT  =  5 

152.  LP 

local  variable 
set  &  used  Sync 

Sync  sets  LP  =  0  and  in  certain  cases  to  LP  =  LP+1 

153.  LPD 

local  variable  pages  (3-66,  3-67) 

set  &  used  Taper(Rgen) 

LPD  =  0,  Chip  counter. 


52 


154.  LRC 


local  variable  pages  (3-94) 

set  &  used  key 

KEY  sets  LRC 

LRC  =  LRC  +  1 

if  LRC  =  NRC  then  LRC  =  0 

LRC  =  1,2,. . .NRC 


155.  LRP 

local  variable  pages  (3-66,  3-67) 

set  &  used  Taper  (Rgen) 

LRP  =  0,  Tape  samples  per  chip  symbol  counter. 


156.  LTAP 

common  block  pages  (3-12  back,  3-17,  3-19,  3-52 

used  detector  to  3-55,  3-59,  3-87,  3-88, 

used  differential  decoder  3-90,  3-95  to  3-98) 

set,  used,  &  output  Sync 
used  Bfilt 

set,  used,  &  output  Parin 
used  Fwate 

Number  of  Backward  Filler  Taps 

LTAP  is  a  fixed  integer  parameter  input  to  the  subroutine  PARIN. 

It  indicates  the  presence  (LTAP>1)  or  absence  (LTAP  =  0)  of  backward 
filter  taos. 

In  the  backward  filter  subroutine,  BFILT,  LTAP  is  sthe  size  of  the 
arrays  JSET  and  BETA  and  is  used  to  define  the  number  of  chip  values 
stored  in  the  array  ALPHA  as  LTAP*NCB. 

The  forward  filter  subroutine,  FWATE,  uses  LTAP  to  define  the  size 
of  the  Kalman  algorithm  arrays  PVAR  and  KVAL  as 

NLTAP=NTAP+LTAP 


157.  LTAP1 

local  variable  pages  (3-85) 

set  &  used  Sync 

LTAP!  =  LTAP/2. 
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158.  LTAPE 

local  variable  pages  (6-1  back,  3-55) 

used  channel 

used  detector 

set  &  used  Taper(Read) 

set  Taper(Rgen) 

used  Sync 

End  of  Input  Tape 

LTAPE  is  the  integer  flag  for  the  end  of  file  reached  on  the  input 
tape.  It  is  initially  set  to  0  by  TAPER-READ  and  set  to  2  when  an  end 
of  file  is  read  on  the  input  tape.  If  TAPER  returns  LTAPE=2  bit  itera¬ 
tion  processing  ends. 

TAPER(RGEN)  sets  LTAPE  =  0 


159.  M 

local  variable  pages  (3-99,  3-100) 

set  Max 

set  &  used  Semul 

M  =  Dimension  of  MSET  (integer,  M  <  N)  where  1  <  M  <  N.  The  final 
value  of  M  depends  on  when  the  ratio  of  the  next  largest  magnitude  to 
the  largest  drops  below  the  dynamic  range  threshold. 


160.  MA 


array 

set  &  used  Max 

The  integer  array 
indicator  array;  MA  is 

MA( I )=0 ,  A( I ) 

MA( I )=1 ,  A( I ) 


pages  (3-99) 

MA(I),  1=1, 2,... N,  is  the  Max- A- Index  Stored 
used  as  follows  in  processing: 

index  I  not  stored  in  MSET 

index  I  stored  in  MSET. 


Initially,  MA(I)=0  for  all  I. 


i 
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161.  MES 


Transmitter  Shift  Register 


array 
set  Pinl 

MES  is  set  by  Pinl  as 
MES(I)  =  MESS(I) 

MES( 1+16)  =  MESS( 1+16) 


162.  MESS 
array 

set,  used,  &  output  Pinl 

MESS  is  set  in  Pinl  to 
MESS(I)  =  M0D(ML ,2) 
MESS( 1+16)  =  MOD  (MH.2) 


163.  MH 

local  variable 
set  &  used  Pinl 

MH  is  set  by  Pinl  as 

MH  =  27945 

MH  =  MH12 


164.  MI 

local  variable  pages  (3-85) 

set  &  used  Sync 

Sync  sets  MI  (delay)  to 

MI  =  | INEXT-IMAXI/NRC  chip  values  and 

MI  =  IMAX- (MSYNC( I )- IMAX ) 
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165.  MODE 


166. 


167. 


168. 


local  variable  pages  (3-56,  3-57) 

used  Parin 

MODE  =  Input  mode  indicator  (Integer). 

Mode  =  0,  Interactive  mode 
Mode  =  1,  Batch  mode. 


MR 

local  variable 
used  Sync 


MSE 

local  variable  pages  (3-48,  3-90,  3-91,  3-96) 

set  DEES 

set,  used,  &  output  detector 
set  &  used  Fwate 
used  Bfilt 

Mean  Square  Error 

DFES  sets  MSE  =  0.0 

MSE  is  the  mean  square  area  real  value  computed  each  bit  symbol 
iteration  for  SM0DE=1  by  the  Detector  as 

MSE  =  1  if  MSE  >  1 

MSE  =  (1.  -  ESTEP)  *  MSE  +  ESTEP  *  TEMPF  *  TEMPF 

Each  iteration  that  the  counter  NUM  is  an  integer  multiple  of  ISKIP  it 
is  used  by  the  Detector  to  compute  ESNR.  It  is  included  as  a  calling 
sequence  argument  to  the  subroutine  BFILT  where  it  is  used  in  the 
decision  skip  the  update  of  the  ALPHA  and  BETA  arrays  if  NALG=2,  KLMS=1, 
and  MSE>0. 10. 


MSELB 

common  block  pages  (5-2  back,  3-12,  3-15,  3-54, 

used  Pinl  3-58) 

set,  used,  &  output  Parin 

Mean  Square  Error  Loop  Bandwidth  in  Hz 

MSELB  is  a  real  input  parameter  to  PARIN  with  a  default  of  1. 

Parin  sets  MSELB  =  1  if  MSELB  =  0 
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169.  MSET 


output  integer  array 


array  pages  (3-99,  3-100) 

set  &  used  Max 

MSET(I)  =  Index  pointer  array  to  A  in  decreasing  magnitude  order, 
1-1,2,. . .M  (Integer) 

A(MSET(1))  =  largest  |A| 

A(MSET(2))  =  2nd  largest  |A| 


A(MSET (M) )  =  Mth  largest  |A| 
integer  array  MSET(l),  1=1,2, ...M 
Max  sets  MSET  to 
MSET(I)  =  J 

where 

J  =  1  ,N 


170.  MSIGN 

local  variable 
set  &  used  Dpgen 

Dpgen  sets  MSIGN  as 

MSIGN  MW0RD(32) 


171.  MSYNC 

array  pages  (3-24  back) 

used  &  output  Sync 

Synchronization  Index  Array 

MSYNC  is  an  integer  array  of  maximum  dimension  MSYNC(180).  It  is 
returned  to  the  subroutine  SYNC  by  MAX  containing  the  index  values  to 
the  synchronization  array  RSYNC  in  order  of  the  decreasing  magnitude 
values  of  RSYNC.  SYNC  uses  MSYNC(l)  and  MSYNC(2)  to  set  IMAX  and  INEXT 
for  the  computation  of  TI  and  NNI. 
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172.  MTAP 


common  block  pages  (3-12,  3-13  back,  3-12,  3-15, 

used  update  input  3- 17,  3-36  to  3-40,  3-55,  3-58, 

used  noise  filter  3-59) 

used  Pinl 

set,  used,  &  output  Parin 
Number  of  Noise  Filter  Taps 

MTAP  is  an  input  parameter  to  the  PARIN  subroutine  with  a  default 
of  1.  MTAP  must  be  an  integer  in  the  range  1<MTAP<25.  If  the  MTAP 
input  is  even  PARIN  forces  it  to  be  odd  by  setting 

MTAP=MTAP+1 

The  value  of  MTAP  then  remains  fixed.  PARIN  also  sets  NFLB=0  for 
MTAP=1. 

Subroutine  Pinl  uses  it  to  set  the  constant  NV  used  in  the  Noise 
Filter  subprogram  as 

NV=(MTAP+l)/2 

MTAP  defines  the  size  of  the  noise  filter  arrays  V(l)  and  XFREQ(K.I), 
1=1 ,2 , . . . MTAP.  When  MTAP=1  the  Noise  Filter  subprogram  is  bypassed. 


173.  MULT 

local  variable 
set  &  used  channel 

MULT  is  set  in  channel  to 

MULT  =  (JBS-D/NTB 
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array 

set  A  used  Dpgen 
Message/Cl  ip  Shift  Register 


pages  (3-69,  3-70) 


MWORD  is  the  shift  register  input  to  the  random  number  generator 
for  the  message  and  chip  data  subroutine  DPGEN. 

Dpgen  sets  MWORD 

MWORD(K)  =  MWORD(K-l) 

MWORD(l)  =  0 

MWORD(K)  =  MWORD(K)  +  JWORD(K) 
if  MWORD(K)  =  2  then  MWORD(K)  =  0 


a)  local  variable 
used  Max 
set  &  used  key 
set  Taper(Read) 


pages  (3-99) 


N  =  Dimension  of  A  (positive  integer) 


b)  array 

used  channel 
set  &  used  Noise 


pages  (3-24  back,  3-22,  3-23,  3-25, 
3-26,  3-79,  3-80) 


Discrete  Noise  Sequence 

N  is  the  complex  noise  sequence  array  returned  by  the  subroutine 
NOISE  dimensioned  N(90).  Indexing  is  N(I),  1-1,2, — NR.  It  is  used  in 
the  Channel  subprogram  in  forming  the  receiver  sample  array  R. 

Noise  sets  N  to 

N( I )  =  N(I-NTB) 

N(I)  =  NDATA(NOP)*SF 

N(I)  =  SUM1*  CMPIX(S(1) , S( 2 ) ) 


176.  NA 


local  variable  pages  (4-5  back,  3-12,  3-14,  3-24, 

used  initialize  parameters  3-53,  3-58,  3-81,  3-82) 

set,  used  &  output  Sync 
used  Pinl 
set  Parin 

Number  of  Bits  Averaged  before  Acquisition  Decision 

NA  is  the  integer  input  parameter  to  initialize  parameters  for  the 
number  of  bit  symbol  iterations  used  by  SYNC  to  reach  an  acquisition 
decision.  For  NA=0  the  Initialize  Parameters  subprogram  sets  SM0DE=1, 
indicating  that  SYNC  is  not  to  be  used.  When  NA>0,  SYNC  decrements  NA 
by  1  each  time  it  is  called  until  NA=0  when  SYNC  computes  NNI ,  TI,  GCON, 
and  sets  SM0DE=1.  Whenever  SM0DE=1  NA  remains  an  unused  constant. 

Parin  sets  NA 
NA  =  0 

if  JTIME  =  1  then  NA  =  20  or  some  chosen  number 
NA  =  Number  of  bit  symbols  to  be  averaged  in  the  correlation  (Integer). 
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177.  NALG 


178. 


co.nmon  block  pages  (3-2  back,  3-12,  3-16,  3-58, 

set  &  used  initialize  parameters  3-87,  3-88) 

used  update  input 

set  channel 

set  &  used  detector 

used  differential  decoder 

used  Fwate 

set  &  used  Parin 

Algorithm  Indicator  (NALG) 

NALG  is  the  integer  indicator  for  the  adaptation  algorithm  used  by 
FWATE  in  updating  the  forward  filter  weight  vector  w.  The  subroutine 
PARIN  sets  NALG  to  the  following  values  according  to  the  input  value  of 
ALGOR: 


ALGOR 

LMS 

KAL 

FIX 

RAKE 


NALG  Weight  Adaptation  Method 

1  Least  Mean  Squares  algorithm 

2  Kalman  algorithm 

3  Weights  remain  fixed  to  initialization  values 

4  Rake  equal i izer 


NALG  remains  constant  unless  the  Kalman  algorithm  is  specified 
(NALG  =  2)  and  the  synchronization  subroutine  is  not  used  (SMODE  =  1). 
When  this  occurs,  the  Initialize  Parameters  subprogram  temporarily  sets 
NALG  =  3  to  force  a  delay  in  the  adaptation.  When  2*NTAP  iterations 
have  been  completed,  NALG  is  reset  to  its  original  input  value  of  2. 


The  Detector  sets  AHAT  =  A  when  NALG  =  2  and  REF  =  3.  The  Detector 
also  calls  the  BFILT  subroutine  to  update  BETA  for  NALG  =  1  or  2  and 
LTAP  >  1. 


Initialize  parameters  sets  NALG  =  3  if  KALA  >  0. 

The  channel  sets  NALG  =  ALGOR  if  KALA  >  0  and  RNUM  >  KALA. 


The  detector  sets  NALG  =  ALGOR. 


NBLOCK 

local  variable 
set  Taper(Read) 

NBLOCK  is  set  by  Taper(Read)  as  NBLOCK  =  5000 
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2^2*. — m. — i. 


local  variable 
set  Taper  (Read) 


pages  (3-63) 


NBYTE  =  3840  Number  of  bytes  per  block  on  the  signal  input  tape. 


NCB 

common  block 
output  update  input 
used  Taper  (Rgen) 
used  Bf i It 
set  Parin 

Number  of  Chip  Symbols  per  Bit  Symbol 

NCB  is  a  fixed  integer  defined  by  Parin  as 
NCB  =  CSR/BSR 

The  BFILT  subroutine  uses  NCB  to  the  ALPHA  array  size  to  LTAP*NCB,  where 
NCB  is  the  number  of  new  ALPHA  values  computed  for  each  bit  symbol 
i terati on. 


pages  (3-16  back,  3-13,  3-19,  3-96, 
3-97) 


NCTAP2 

local  variable 
set  &  used  Parin 

Parin  sets  NCTAP2  to 

NCTAP2  =  2*CTAP 


NFLB 

common  block  pages  (3-7  back,  3-12,  3-18,  3-36, 

used  initialize  parameters  3-37,  3-39,  3-40,  3-55,  3-58) 

set  &  used  noise  filter 

used  detector 

set,  used,  &  output  Parin 

Noise  Filter  Loop  Bandwidth  in  Hz 

NFLB  is  a  real  input  parameter  to  the  subroutine  PARIN.  The 
Initialize  Parameters  subprogram  uses  it  in  computing  the  number  of 
iterations  for  NOISE  FILTER  adaptation  IX  and  the  Noise  Filter  sub¬ 
program  uses  it  in  computing  VSTEP. 

When  NFLB=0,  the  noise  filter  weights,  V(I),  I=1,2,--MTAP,  are 
fixed.  If  NFLB>0,  NFLB  is  reset  to  0  when  IBS  reaches  IX.  The  sub¬ 
routine  PARIN  sets  NFLB=0  when  MTAP=1. 


183.  NGEN 


local  variable  pages  (4-6  back,  3-62,  3-63,  3-65, 

set  &  used  update  input  3-66) 

used  Taper  (Read) 
used  Taper  (Rgen) 

Number  of  Input  Samples  Accessed  by  TAPER 

NGEN  is  the  integer  number  of  new  input  samples  for  TAPER-RGEN  to 
generate  or  the  number  of  tape  samples  for  TAPER-READ  to  generate.  It 
is  a  calling  sequence  argument  for  both  versions  of  TAPER.  NGEN=NTB+NRQ 
for  each  bit  symbol  iteration  unless  it  is  the  first  iteration  where 
SM0DE=1.  (KSYNC=1  and  SM0DE=1).  Then  it  is  increased  by  NSHIFT. 


184.  NI 

local  variable 
set  &  used  Sync 

NI  =  1-1  as  set  by  Sync 


185.  NIN 

Noise  input  indicator 

local  variable 
set  DFES 

set  &  used  initialize  parameters 
used  differential  decoder 
used  Noise 
set  &  used  Parin 

NIN  =  0  as  set  by  DFES,  Initialize  Parameters,  and  Parin 
Parin  also  sets  NIN  =  1 


186.  NITAP 

local  variable  pages  (3-90,  3-91) 

set  &  used  Fwate 

is  NlTAP  =  NTAP  +  LTAP  where  NTAP  is  the  number  of  forward  filter  taps 
and  LTAP  is  the  number  of  backward  filter  taps. 
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local  variable 
set  &  used  Noise 


NN  is  set  by  Noise  to 
NN  =  NTB 

NN  =  NR  if  RIBS  =  1 


NNI 

common  block  pages  (3-9  back,  3-11,  3-17,  3-19, 

set  DFES  3-56,  3-58,  3-63,  65,  3-66, 

output  update  input  3-83,  3-85) 

set  channel 

used  Taper  (Rgen) 

set  &  output  Sync 

set  Parin 

Number  of  Input  Samples  to  Be  Skipped  in  Order  to  Produce  Bit 
Synchroni zati on 

When  the  SYNC  subroutine  is  required  for  acquisition  (NA  >  0  and 
SMODE  =  0),  NNI  is  determined  by  SYNC.  The  initial  value  of  NNI  is  0. 

If  SYNC  is  not  required  (NA  =  0),  NNI  is  an  input  parameter  to  PARIN. 

The  value  of  NNI  set  by  SYNC  or  input  by  PARIN  is  used  for  only  one  bit 
symbol  iteration  call  (NUM  =  1)  to  the  Tape  Read  subroutine  TAPER  (Rgen) 
After  each  call  to  TAPER,  NNI  is  reset  to  0. 

TAPER  (Rgen)  uses  NNI  in  computing  the  number  of  input  samples  to 
read  as 


NK  =  NGEN  +  NNI 


NNI  is  set  to  NNI  =  0  by  Sync,  DFES,  and  channel  Sync  sets 


NNI  =  I MAX-1 


189.  NOSIG 


TT3T 


common  block  pages  (3-18  back,  3-17,  3-22,  3-26, 

used  and  output  initialize  parameters  3-55,  3-57) 

used  channel 

set  &  used  noise  filter 

used  differential  decoder 

Set  &  output  Parin 

l  Input  Signal  Control  Indicator 

NOSIG  is  an  input  integer  parameter  of  0  or  1  to  PARIN.  N0SIG=0 
is  the  normal  operation  mode.  When  NOSIG  =  1,  there  is  no  signal  input, 
.  R(K)  =  N(K).  If  NFLB  =  0,  PARIN  sets  NOSIG  =  0. 

The  N0SIG-1  option  is  no  longer  used  in  the  program. 


190.  NOZDC 

common  block  pages  (3-21  back,  3*22,  3-26, 

used  channel  3-55,  3-57) 

set  &  output  Parin 

DC  Noise  Indicator 

NOZDC  is  a  fixed  integer  parameter  input  to  PARIN.  N0ZDC=0  is 
the  normal  default  case.  NOZDC-1  implies  a  test  case  of  DC  NOISE. 


191.  NPD 

local  variable  pages  (3-66) 

set  &  used  Taper  (Rgen) 

NPD  =  CSR/BSR,  Number  of  chip  symbols  per  bit  symbol 
NPD  =  NCB 


192.  NPN 

local  variable  pages  (3-81,  3-82) 

set,  used,  &  output  Sync 

NPN  =  RANGE*6000*R$R*1. E-9  =  QN 

NPN  =  NPN+1  IF  QN  >  NPN 

The  program  exits  if  NPN  is  greater  than  180 


193.  NPOW 


local  variable 
set  &  used  channel 
used  &  output  detector 

NPOW  =  JBS  -  MULT  *  NTB 


194.  NQ 

local  variable  pages  (3-15) 

set  initialize  parameters 
output  update  input 
used  channel 

NQ  =  NTB+INT-NTR+KSET(CTAP) 


195.  NR 

common  block 
output  update  input 
used  Taper  (Both) 
used  Noise 

set,  used,  &  output  Parin 
used  channel 

Number  of  Receiver  Input  Samples  -  Integer 

NR  is  the  size  of  the  Receiver  input  array  R. 

Parin  sets  the  fixed  value  as 

NR  =  NTB+INT-NTR 


196.  NRAN1 

local  variable 
set  DFES 
output  channel 
used  Gauran 
set  Parin 

NRAN1  is  set  by  DFES  and  Parin  as 
NRAN1  =  0 


pages  (3-9  back,  3-13,  3-19,  3-63, 
3-80) 
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NRAN2 

local  variable 
set  OFES 
output  channel 
used  Gauran 
set  Parin 

NRAN2  is  set  by  DFES  and  Parin  as 
NRAN2  =  0 


NRB 

common  block 

used  initialize  parameters 
used  &  output  update  input 
used  interpolator 
used  noise  filter 
used  key 
used  sync 
used  Fwate 
set  &  used  Parin 

Number  of  Receiver  Samples  per  Bit  Symbol 

NRB  is  set  to  the  fixed  integer  value 

NRB=RSR/BSR 

It  is  used  by  the  Initialize  Parameters  Subprogram  to  compute  the 
size  NY  of  the  forward  filter  output  array  Y  and  by  the  Update  Input 
subprogram  to  compute  the  size  NS  of  the  interpolator  input  array  S 

NY=NRB 

NS=NRB=NSPAN 

In  the  noise  filter  subprogram  it  is  the  number  of  XFREQ  (K,I), 

K  =  1,2, — NRB,  values  generated. 

The  SYNC  subroutine  uses  NRB  to  set  the  size  of  the  ALPHA  array 

LIMIT=NPN+NRB-1 

where  NRB  is  the  number  of  new  ALPHA  values  generated. 


pages  3-11  back,  3-13,  3-19,  3-37 
to  3-39,  3-41  to  3-43,  3-81, 
3-82,  3-87,  3-89,  3-93,  3-94) 


The  number  of  PN  sequence  values  to  be  generated  by  KEY  is  I  =  NRB 
and  the  number  of  interpolator  output  x  values  used  in  computing  RSYNC 


199.  NRC 


common  block  pages  (3-17  back,  3-13,  3-19, 

output  update  input  3-93,  3-94,  3-97) 

used  Sync 

set  Parin 

used  Bfilt 

used  Key 

Number  of  Receiver  Samples  per  Chip  Symbol 

The  constant  NRC  is  defined  in  Parin  as  NRC=RSR/CSR.  It  is  used 
by  the  KEY  subroutine  as  the  number  of  PN  sequence  duplications. 

The  BFILT  subroutine  uses  NRC  as  part  of  the  PN  index  in  computing 
the  ALPHA  and  GBACK  arrays. 


200.  NRD 

local  variable  pages  (3-63,  3-66) 

set  &  used  Taper  (Both) 

For  IBS=1  only,  the  number  of  forced  delay  samples  is  set  (NRD=5). 
Later  set  to  NRD=NRD-1 


201.  NREC 

local  variable 
used  DFES 


202.  NRP 

local  variable  pages  (3-66) 

set  &  used  Taper  (Rgen) 

NRP  =  NTR*RSR/C5R,  Number  of  tape  samples  per  chip  symbol. 


203.  NRQ 

local  variable 

set  &  used  update  input 

used  chanel 

Update  input  sets  NRQ  as 

NRQ  =  0 
and 

NRQ  =  INT  -  NTR 
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common  block  pages  (3-14  back,  3-16,  3-19) 

set  DFES 

set,  used,  &  output  update  input 
used  interpolator 
used  Fwate 

Interpolator  Array  Size 

NS  is  the  fixed  integer  dimension  of  the  Interpolator  arrays  S  and 
X.  The  Update  Input  subprogram  sets 

NS=NRB+NSPAN 

NX=NS 

NS  has  the  range  1  <  NS  <  500  as  set  in  DFES 


NSHIFT 

local  variable  pages  (3-18,  3-42) 

set  DFES 

set,  used,  &  output  update  input 

The  NSHIFT  parameter  is  computed.  The  purpose  of  this  parameter  is 
to  maintain  synchronization  for  the  nondispersive  channel  for  any 
forward  filter  tap  specification. 

NSHIFT  =  NTR*ISET ( 1) . 

NSHIFT  is  set  to  0  by  DFES. 


NSPAN 

common  block  pages  (3-9,  3-10  back,  3-18,  3-19) 

set  DFES 

set,  used,  &  output  update  input 
Forward  Filter  Span 

In  receiver  sample  widths  DFES  sets  NSPAN  =  0 
NSPAN  is  an  integer  set  to  the  fixed  value 
NSPAN  =  MJX  I SET  (I)  =  ISMAX 

by  the  update  Input  subprogram  where  it  is  used  to  compute  the  size  NS 
of  the  Interpolator  arrays  S  and  X 


NS=NRB+NSPAN 


207.  NTAP 


common  block 

used  initialize  parameters 
used  update  input 
used  forward  filter 
used  differential  decoder 
set,  used  &  output  Sync 
set,  used,  &  output  Parin 
used  Fwate 
used  Bf i 1 t 

Number  of  Forward  Filter  Taps 

NTAP  is  a  fixed  integer  parameter  to  the  subroutine  PARIN  with  a 
default  of  1.  It  defines  the  size  of  the  forward  filter  arrays  ISET  and 
W.  The  subroutine  FWATE  uses  NTAP  to  set  the  size  of  the  Kalman 
algorithm  arrays  VKAL  and  PVAR. 

NLTAP=NTAP+LTAP 

FWATE  also  uses  NTAP  as  an  index  indicator  in  computing  the 
variable  CON  for  NALG=2. 

The  Backward  Filter  subroutine  uses  NTAP  in  the  index  for  the  VKAL 
array  when  KLMS=0. 

The  Initialize  Parameters  subprogram  uses  NTAP  when  NALG=2  and 
SM0DE=1  to  set  the  adaptation  delay  to  KALA=2*NTAP  iterations. 

Sync  sets  NTAP  =  NTAP  +  1  i f  MOD  (NTAP,  2)  *  0. 


pages  (3-13 

back,  3-13,  3-16, 

3-20, 

3-41, 

3-52,  3-54,  3-55, 

3-58, 

3-59, 

3-83,  3-84,  3-87, 

3-89, 

3-90, 

3-92,  3-96,  3-98) 

208.  NTAP1 

local  variable 
set  &  used  Sync 

Sync  sets  NTAP1  to  be 

NTAP1  =  NTAP/2 


209.  NTAP2 

local  variable 
set  &  used  Parin 

Parin  sets  NTAP2  to  be 

NTAP2  =  2  *  NTAP 
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210.  NTB 

Number  of  channel  samples  per  bit  symbol 

common  block  pages  (3-10,  3-11  back,  3-13,  3-15, 

used  initialize  parameters  3-19,  3-23,  3-24,  3-80) 

used  &  output  update  input 

used  Noise 

set  &  used  Parin 

used  channel 

Number  of  Channel  Samples  per  Bit  Symbol 

NTB  is  a  fixed  integer  value  set  by  Parin  to 
NTB=TSR/BSR 

where  it  is  used  to  set 
NR=NTB+INT-NTR 

and 

NQ=NTB+INT-NTR+KSET(CTAP) 

Each  bit  symbol  iteration  is  used  to  set  the  number  of  receiver  samples 
NGEN  to  be  input  by  TAPER.  The  Channel  subprogram  uses  NTB  in  computing 
the  index  NPOW  and  the  NOISE  subroutine  uses  it  to  set  the  number  of 
noise  samples  generated  in  N  array. 


211.  NTR 

common  block  pages  (3-13,  3-15,  3-18,  3-42,  3-80) 

used  initialize  parameters 

used  Taper  (Rgen) 

used  &  output  update  input 

used  interpolator 

set  &  used  Parin 

NTR  =  TSR/RSR  as  set  by  Parin 
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212.  NV 


local  variable  pages  (3-15,  3-63  to  3-65) 

used  noise  filter 

set  &  used  Taper  (Read) 

set  Pin  1 

Taper  (Read)  sets  NV  to 
NV=NG£N+NNI 

The  number  of  data  samples  to  transfer  to  the  output  array,  VR, 
plus  the  number  of  data  samples  to  be  skipped  for  bit  synchroni¬ 
zation  and  NV  =  NV  -  1. 

An  integer  half-width  of  the  Noise  Filter  is  computed 
NV  =  (MTAP+l)/2  as  set  in  Pin  1. 


213.  NVIN 

common  block  pages  (3-13  back,  3-63,  3*64) 

set  &  used  Taper  (Read) 

Receiver  Sample  Pointer 

NVIN  is  an  integer  variable  used  by  the  Tape  Read  subroutine  TAPER 
(file  name  READ)  as  the  pointer  to  the  consecutive  receiver  samples  in 
the  input  data  block  array  VIN.  When  IBS=1,  NVIN  is  initialized  to  0 
in  Taper  (Read). 

The  data  sample  pointer  is  set  (NVIN=0). 

Later  updated  NVIN  =  NVIN  +  1 


214.  NX 

common  block  pages  (3-12  back,  3-19,  3-37,  3-39 

set  DFES  to  3-42,  3-87) 

set  update  input 

used  interpolator 

used  noise  filter 

used  forward  filter 

Noise  Filter  Output  Size 

NX  is  the  fixed  integer  size  of  the  noise  filter  output  array  x. 
The  Update  Input  subprogram  sets 


NX=NS 


OFES  sets  NX  initially. 
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215.  NY 

local  variable 
set  initialize  parameters 
used  forward  filter 
used  Key 
used  compressor 

Forward  Filter  Output  Size 

The  integer  size  NY  of  the  forward  filter  output  array  Y  is  set  by 
the  Initialize  Parameters  subprogram 

NY=NRB . 

It  is  used  by  the  Compressor  as  a  calling  sequence  argument  for  KEY 
to  generate  NY  PN  sequence  values. 


216.  ONE 

local  variable 
set  &  used  ERFC 

ONE  is  set  by  ERFC  as 

1.D0  in  a  Data  Statement 


217.  P 

local  variable  pages  (3-67,  3-93,  3-94) 

set  &  used  ERFC 

set  &  used  key 

set  &  used  Taper  (Rgen) 

The  receiver  sample  P  (Complex)  is 

P  =  P*EJ. 

ERFC  sets  the  value  of  P  in  a  Data  Statement. 


pages  (5-3  back,  3-13,  3-44,  3-93, 
3-94) 
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218.  PI 


common  block 

used  initialize  parameters 
set  DFES 
used  channel 
used  noise  filter 
used  Bfilt 

n 

The  real  constant  PI  is  set  by  the  DFES  main  program  to 
PI  =  3. 14159265. 


pages  (3-7  back,  3-14,  3-36,  3-92 
3-109) 


219.  P1X 

local  variable 
set  &  used  Sine 

PIX  is  set  by  Sine  as 

PIX  =  PI  *  X 


220.  PN 

array  pages  (3-24  back,  3-36,  3-39,  3-81, 

used  noise  filter  3-82,  3-88,  3-89,  3-94,  3-96, 

used  Sync  3-97) 

used  Fwate 

used  Bfilt 

set  Key 

Pseudo-Noise  Sequence  Array 

The  complex  array  PN(K),  K-1,2,--NY  generated  each  bit  symbol 
iteration  by  the  subroutine  KEY.  It  is  a  calling  sequence  argument  of 
the  subroutines  KEY,  SYNC,  BFILT,  and  FWATE. 

The  Noise  Filter  subprogram  uses  PN  in  computing  the  array  V.  The 
Compressor  subprogram  uses  it  to  compute  the  compressor  output  Z.  SYNC 

uses  PN  to  compute  the  array  ALPHA  and  FWATE  to  compute  the  S  BAR  array 

PN  is  used  in  BFILT  for  setting  the  ALPHA  and  GBACK  arrays. 


Key  sets  PN(N)  =  CONJG(P) 


221. 


222. 


223. 


POW 

array  pages  (3-25  back,  3-13,  3-17,  3-22 

set  DFES  >54,  3-60) 

used  channel 

set  &  output  Parin 

Channel  Relative  Power  Array 

POW  is  a  real  array  of  the  subroutine  PARIN  dimensioned  P0W(4). 

Its  default  is  P0W(1)=1,  P0W(I)=0,  I=2,3,CTAP.  It  is  the  relative  power 
of  path  I  with  respect  to  the  first  path.  It  is  used  by  the  Channel 
subprogram  for  the  variance  CVAR  used  in  the  channel  random  number 
generation  by  GAURAN 

CVAR  =  P0W( J )/2 


PSK 

common  block  pages  (3-2  back,  3-45,  3-46,  3-48, 

used  detector  3-54,  3-57) 

set,  used,  &  output  Parin 

Number  of  Transmitted  Phases 

PSK  is  an  integer  of  fixed  value  2  or  4.  It  is  used  by  the 
Detector  as  an  indicator  to  determine  the  value  of  AHAT.  PSK  is  an 
input  parameter  to  the  subroutine  PARIN  with  the  default  PSK  =  2. 


PVAR 

virtual  array  pages  (3-89  to  3-91) 

set,  used,  &  output  Fwate 

The  inverse  matrix  estimate  PVAR 

The  Kalman  matrix  PVAR  is  updated. 

PVAR  is  set  by  Fwate  to 

PVAR(I,J)  =  (4. ,0. ) 

if  I  =  J  PVAR(I,J)  =  (1. ,0. ) 

PVAR(I,I)  =  PVAR( I , I )*AKC-RCON*(CABS(VKAL( I ) )**2) 

PVAR(I,J)  =  PVAR( I , J)*AKC-RC0N*VKAL( I )*C0NJG(VKAL(J)) 

PVAR( J , I )  =  C0NJG(PVAR(I,J)) 
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224.  Q 


virtual  array  pages  (3-22  to  3-24,  3-27) 

set  &  used  channel 

Every  bit  symbol  iteration,  the  subprogram  receives  a  transmit 
signal  array 

Q(K)  =  Input  array  to  the  subprogram  (Complex). 

K  =  1 ,  2 ,  ...  NQ .  NQ  <  330. 
NQ=NTB+INT-NTR+KSET(CTAP) . 

the  most  recent  input  values  (Q  array) 

Q(K)  =  Q(K-NTB)  K  =  NTB  +  1,  ...  NQ 

Q(K)  =  CMPLX(VR( 1,K),VR(2,K) 

225.  QE 

local  variable 
set  &  used  detector 

QE  is  set  by  the  detector  as 

QE  =  AIMAG(E) 

QE  =  SIGN  (1. ,  QE) 


226.  QN 

local  variable 
set  &  used  Sync 
set  Parin 


QN  is  set  by  Sync  as 

QN  =  6000.  *  RANGE  *  (l.E-9)  *  RSR 
QN  is  set  by  Parin  as 

QN  =  6000.  *  (1.0  E-9)  *  RSR 
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227.  QS 

local  variable 

set  &  used  interpolator 

The  interpolator  secs  QS  as 

QS  =  (0.0,  0.0) 

QS  =  QS  +  GCON  *  R  (NTR  *  K  -*■  I  +  INTI  +  ( 1-NTR) )  *  SINK(IQ) 

228.  QW 

local  variable 
set  &  used  Fwate 
set  &  used  Bfilt 

QW  is  set  by  Fwate  as 

QW  =  AIMAG(W( I ) ) 

QW  =  100.  *  SIGN(l.QW) 

QW  is  set  by  Bfilt  as 

QW  =  A I MAG  (BETA(I)) 

QW  =  100.  *  SIGN( 1. ,  QW) 

229.  R 

array 
set  OFES 

set,  used,  &  output  channel 
used  interpolator 
output  detector 

R  is  the  receiver  input  array 

DFES  sets  (R(I)  =  (0.0,  0.0) 

(Channel  sets  R  to 

R(K)  =  (0.0.) 

R(K)  =  R(K)  +  H( I )  *  Q  (K-VSET(I)) 

R(K)  =  R(K)  +  N(K) 
where  K  =  1,  NTB  +  NRQ 


pages  (3-14,  3-23  to  3-27,  3-29, 
3-30  to  3-33) 
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230.  RANGE 


common  block  pages  (3-7  back,  3-55,  3-58,  3-81, 

set  DFES  3-86) 

used  &  output  Sync 

set,  used,  &  output  Parin 

Range  in  Nautical  Miles 

RANGE  is  an  input  parameter  to  PARIN  only  if  the  bit  synchroniza¬ 
tion  routine  SYNC  is  to  be  used  (NA  >  0).  It  is  used  by  SYNC  to  compute 
the  number  of  PN  sequence  values  range  used  in  synchronization  NPN.  The 
default  is  RANGE  =  300.  DFES  initially  sets  RANGE  =  0.0. 


231.  RATE 

local  variable  pages  (3-52) 

set,  used,  &  output  differential  decoder 

RATE  -  Bit  error  rate  (Real). 

The  error  rate  (RATE)  is  computed  as  the  number  of  bit  errors 
divided  by  the  total  number  of  received  bits. 

The  differential  decoder  sets 

RATE  =  ERR0R/(RNUM  *2) 

and  if  REF  =  1  then 

RATE  =  RATE  *  2 


232.  RCON 

local  variable  pages  (3-91) 

set  &  used  Fwate 

RC0N=KSTEP*AKC . 


233.  RDEL 

local  variable 
used  detector 
set  Pinl 

Pin  1  sets  RDEL  as 

RDEL  =  2.  *  PI  *  RGLB/BSR 
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nmtaa 


local  variable 
set  &  used  detector 

RE  is  set  in  the  detector  to 

RE  =  REAL(E) 

RE  =  SIGN( 1. ,  RE) 


common  block 

set  update  input 

used  detector 

used  differential  decoder 

set,  used,  &  output  Parin 

set  Fwate 


pages  (3-2  back,  3-11,  3-15,  3-45  to 
3-50,  3-54,  3-57,  3-92) 


Presence  or  Absence  of  Reference  Signal  Indicator 

REF  is  an  integer  input  parameter  to  PARIN  of  value  0,  1,  2,  or  3 

PARIN  uses  the  input  value  of  REF  to  set  the  differential  decoder 
indicator  IDEC.  The  value  of  REF  remains  constant  for  NALG  =  1  or 
NALG  =  3.  When  NALG  =  2,  the  subroutine  FWATE  modifies  REF  in  the 
Kalman  algorithm  adaptation  to  the  following: 

REF  =  3  for  KNUM  =  KVAR 

REF  =  0  for  KNUM  =  KENO 

When  REF  =  3  and  NALG  =  2,  the  Detector  sets  AHAT  =  A.  The  Detector 
also  uses  the  current  value  of  REF  as  an  indicator  for  setting  the 
adaptation  error  sample  E. 


RGAIN 

local  variable 
set  DFES 

set,  used,  &  output  detector 
RGAIN  is  set  by  DFES  as 
RGAIN  =  (1.0,  0.0) 

RGAIN  is  set  in  the  detector  to 

RGAIN  =  RGAIN  +  RDEL  *  E  *  CONJG(Z) 


237.  RGLB 

Rake  gain  loop  bandwidtn  in  HZ,  default  =  5 
common  block 

set,  used,  &  output  Parin 
used  Pinl 

Parin  sets  RGLB  to 

RBLB  =  5.0  or  some  other  number 
if  RGLB  =  0  then  RGLB  =5.0 


238.  RIBS 

Number  of  bit  Symbol  iterations  including  Sync 

common  block 
set  DFES 

set,  used,  &  output  update  channel 

used  interpolator 

used  detector 

used  noi se  f i 1  ter 

used  Noise 

used  Key 

set,  used,  &  output  Sync 
set  Taper  (Read) 
used  Taper  (Rgen) 

RIBS  is  set  as 

RIBS  =  0.0  in  DFES 

Update  input  sets  RIBS  to 

RIBS  =  RIBS  +  1 


239.  RJBS 

local  variable 
set  &  used  channel 

Channel  sets  RJBS  to 

RJBS=RIBS 

If  RJBS  >  l.E+4  then  RJBS  =  1 
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240.  RKAL 


241. 


242. 


243. 


244. 


local  variable 
used  update  input 
set  Pinl 

RKAL  is  set  by  Pin  1  as 

RKAL  =  FLOAT  (KRST/KADAPT) 


RKBS 

local  variable 

set  &  used  noise  filter 

The  noise  filter  sets  RKBS  as 

RKBS  =  RIBS 

RKBS  =  RNUM  if  SMODE  =  1 


RKIBS 

local  variable 
used  forward  filter 
set  &  used  detector 


RKIPS 

local  variable 
set  DFES 

RKIPS  is  set  by  DFES  to 
RKIPS  =  0.0 

RMAX 

local  variable 

set,  used,  &  output  Sync 

Sync  sets  RMAX  to 

RMAX  =  0. 

RMAX  =  RS(IMAX) 
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245 . 


246. 


247. 


248. 


1 


RNBS 

Real  number  of  bit  symbol  iterations  to  do  after  Sync 
common  block 

used  differential  decoder 
set  &  output  Parin 

RNBS  is  set  by  Parin  to  a  number  desired  for  the  number  of  bits, 
example: 


1,000,  10,000,  etc. 


RNEXT 

local  variable 
set  &  output  Sync 

Sync  sets  RNEXT  to 

RNEXT  =  0 

RNEXT  =  RS(INEXT) 


RNUM 

Number  of  bit  symbol  iterations  after  Sync 

common  block 
set  DFES 

set  update  input 
set  &  used  channel 
used  noise  filter 
used  detector 

used  &  output  differential  decoder 
used  Bfilt 

RNUM  is  set  RNUM  =0.0  in  DFES 

and 

Update  Input  Channel  sets  RNUM  to 
RNUM  =  RNUM  +  1 


RS 

array  pages  (3-83) 

set,  used,  &  output  Sync 

Comparision  is  made  on  the  magnitude  of  the  correlator  output 
RS( I )  =  CABS(RSYNC(I)) 
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249.  RSR 


common  block 

used  initialize  parameters 
used  noise  filter 
used  Taper  (Rgen) 
used  Sync 

set,  used,  &  output  Parin 
used  Pin  1 

Receiver  Sample  Rate  in  Hz 

RSR  is  a  real  input  constant  to  the  PARIN  subroutine.  It  has  a 
default  value  of  96K  and  must  be  an  integer  multiple  of  BSR. 

RSR  is  used  by  the  Initialize  Parameters  subprogram  in  computing 
NRB,  NRC,  NTR,  and  the  noise  filter  adaptation  delay  IX. 

In  the  Noise  Filter  subprogram  it  is  used  in  computing  the  constant 


pages  (3-4  back,  3-12,  3-24,  3-29, 
3-30,  3-36  to  3-38,  3-51, 
3-55,  3-58,  3-66,  3-72,  3-80, 
3-86) 


ALPHA  =  J  (lfl) 


and  the  variable 


Cx  =  2nd  ‘  I  =  1,  2,  -  MTAP 


When  all  bit  symbol  iterations  have  been  completed,  RSR  is  used  in 
computing  the  measured  SNR 

SNRM  =  10  log1Q  (SP0W*TSR/(*BSR*(SIGMA**2))) 

The  SYNC  subroutine  uses  RSR  in  computing  the  number  of  PN  sequence 
range  values  used  in  bit  synchronization 

NPN  =  6000*RANGE*RSR*10"9 
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250.  RSYNC 


251. 


252. 


array  pages  (3-25  back,  3-82,  3-83) 

set  it  used  Sync 

Bit  Synchronization  Array 

RSYNC  is  a  complex  array  computed  by  SYNC  of  size 

NPN=6000*RANGE*10_9*RSR  where  NPN  <  180. 

SYNC  uses  it  as  the  calling  sequence  input  array  for  the  subroutine 

MAX. 

SYNC  sets  RSYNC(I)  =  (0.  ,0.  ) 

RSYNC(I)  =  RSYNC(l)  +  X(J)  *  ALPHA(J+NI) 

RTHQLD 

common  block 

used  forward  filter 

set,  used,  &  output  Parin 

Rate  threshold  in  dB  down  from  maximum  weight,  default  =  12 
Parin  sets  RTHOLD  as 

RTHOLD  =  10 . **(-RTHOLD/20. ) 

and 

if  RTHOLD  =  0  then  RTHOLD  =  12 


RW 

local  variable 
set  &  used  Fwate 
set  &  used  Bfilt 

RW1S  set  by  Fwate  as 

RW  =  REAL(W(I)) 

RW  =  100. *SIGN(1. , RW) 

RW  is  set  by  Bfilt  as 

RW  =  REAL  (BETA(I)) 

RW  =  100. *SIGN(1. , RW) 
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253.  S 

Interpolator  output  array 

a)  Virtual  array  pages  (3-16,  3-30  to  3-33,  3-37) 

set,  used,  &  output  interpolator 

used  noise  filter 

The  interpolator  sets  S  as 

S(K)  =  S(K-NRB)  where  K  =  NS,NRB+1,-1 
S(K)  =  QS  where  K  =  1,  NRB 

b)  array  pages  (3-105) 

set  Gauran 

used  Noise 

The  real  and  imaginary  parts  of  the  complex  number  are  then 
formed 

S(l)  =  (-2  V  loge  A)*5  cos  (2ttB) 

S ( 2 )  =  (-2  V  loge  A)5*  sin  (2nB) 

S(l)  =  Real  part  of  the  complex  number 
S(l)  =  Imaginary  part  of  the  complex  number 
Gauran  sets  S  to  S  corresponds  to  XC 

S(l)  =  X*C0S(Y) 

S(2)  =  X*SIN(Y) 

254.  SBAR 

array  pages  (3-89) 

set  &  used  Fwate 

weight  input  voltage  after  compression 

NRB  _ 

SBAR(I)  =  I  PN(K)*X(ISET(I)+K) 

K=1 
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255.  SCALE 


256. 


257. 


local  variable 

set  &  used  Taper  (Read) 

Taper  (Read)  sets  SCALE  as 

SCALE  =  1.0E-4 


SDC 

local  variable 
set  noise  filter 

SDC  is  set  in  the  noise  filter  as 
SDC  =  0. 


SI 

local  variable 
set  &  used  detector 

The  detector  sets  SI  as 

SI  =  SIGN(1. 0,  AIMAG(ZC) 
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258.  SIGMA 


local  variable  pages  (5-3,  5-4  back,  3-18, 

set  initialize  parameters  3-49,  3-51,  3-77,  3-79 

used  channel 
set  Noise 

Standard  Deviation  o  of  a  Quadrature  Component  of  the  Complex 
Sample  N 

SIGMA  is  a  real  constant  set  to 
SIGMA  =  SQRT(TSR/(BSR*2. *SNR)) 

by  the  Initialize  Parameters  subprogram  and  NOISE  subroutine 
In  Noise  it  is  set  to 
SIGMA  =  SQRT(SVAR) 

The  Channel  subprogram  uses  it  to  form  the  noise  array  N 
N(K)  =  0+ jo 
when  NOZDC  =  1. 


259.  SINC 

local  variable  pages  (3-109) 

used  Sine 

set  update  input 

The  function  SINC(X)  is  defined  as 

_  sinrtx  ,  ,  n 

SINC(X)  nx  ’  for  x  *  0 

=  1  ,  for  X  =  0. 

For  X  =  0  or  |sinnx|  -  |nx|  <  10  ^  the  routine  returns 
SINC  =  1 
otherwi se 

SINC  =  (SIN(PIX))/ (PIX) 


-25, 

3-80) 


Noi  se 
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260.  SINK 


array  pages  (3-19  3-29,  3-33) 

set  &  output  update  input 
used  interpolator 

SINK  (I)  Array  of  sin  x/x  interpolation  values  (Real). 

I  =  1,2,. . . INT. 


SINK(J)  =  SINC 


. J- INTI- 1 
(  NTR 


-  TI) 


where  the  SINC(X)  function  is  defined  in  (3.4.22.). 

To  speed  up  the  computation,  the  SINC  (•)  function  is  precomputed 
as  an  array  SINC  in  the  subprogram  UPDATE  INPUT  (3.4.2)  as  follows 

SINC(J)  =  SINC( I/NTR-TI ) ,  I  =  -INTI,  ...  INTI 

J  =  I+INT1+1. 


SINK(J)  =  SINC(XX) 


261.  SIX 

local  variable 

set  &  used  initialize  parameters 

SIX  is  set  in  initialize  parameters  to 
SIX  =  3. *RSR/(2. *PI*NFLB*NRB) 


262.  SKIP 

common  block 

used  detector 

used  differential  decoder 

used  channel 

used  noise  filter 

used  &  output  Fwate 

set  &  used  Parin 

SKIP  is  set  in  Parin  to  the  number  of  bits  to  be  skipped  during 
printout  according  to  desired  number,  example: 

SKIP  =  100,  1,000,  or  5,000 
If  SKIP  =  0,  SKIP  =  100 
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263.  SMODE 


local  variable  pages  (6-1  back,  3-12,  3-14,  3-16, 

set  &  used  initialize  parameters  3-81,  3-83,  3-85) 

used  update  input 

used  interpolator 

jsed  channel 

used  noise  filter 

set  Sync 

Bit  Synchronization  Mode 

SMODE  is  the  integer  synchronization  mode  indicator.  Its  values 
are  0  or  1.  When  SMODE  =  0  bit  synchronization  is  performed  by  SYNC  for 
IBS=NA  iterations  until  the  acquisition  decision  is  obtained.  When 
SM0DE=1,  the  acquisition  decision  has  been  made  and  SYNC  is  not  used. 

The  Initialize  Parameters  subprogram  sets  SM0DE=0  for  NA>0  and  SM0DE=1 
for  NA=0.  The  subroutine  sets  SM0DE=1  after  NA  iterations  when  it  makes 
the  acquisition  decision.  SMODE  is  a  calling  sequence  argument  of  the 
subroutine  SYNC. 


264.  SNR 

common  block  pages  (4-6  back,  3*18,  3-53,  3-58, 

set  &  used  initialize  parameters  3-72,  3-74  to  3-78) 

used  Noise 

set  &  output  Parin 

Signal  to  Noise  Ratio 

SNR  is  the  real  variable  input  parameter  of  the  signal  to  noise 
ratio  in  dB  for  the  subroutine  PARIN.  The  Initialize  Parameters  sub¬ 
program  converts  it  to 

SNR=10. **(SNR/10. ). 

SNR  is  a  calling  sequence  argume  t  for  the  NOISE  subroutine  where  it  is 
used  to  compute  the  variance  SVAR  for  the  random  number  generator  GAURAN. 


265.  SNRM 

local  variable  pages  (3-52) 

set  &  output  differential  decoder 

SNRM  =  Measured  signal-to-noise  ratio  in  dB  (Real). 

It  is  also  printed.  The  measured  signal-to-noise  ratio  expressed  in  dB 
is  calculated  as 

SNRM  =  10. *ALOG10(SPAW*TSR/(2. *BSR*(SIGMA**2))) 
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loca1  variable 

set  DFES 

set  &  used  channel 

used  differential  decoder 

DFES  sets  SPAW  to 

SPAW  -  1.0 

The  channel  sets  SPAW  to 

SPAW=( 1. - 1. /RNUM)*SPAW+( (CABS( R(NPOW) ) )**2)/ (2. *RNUM) 


SPOW 

local  variable  pages  (3-14,  3-24,  3-25,  3-49,  3-51) 

set  DFES 

set  &  used  channel 

output  detector 

used  differential  decoder 

DFES  sets  SPOW  as 

SPOW  =  1. 

Receiver  input  power  per  quadrature  channel  and  recursively 
computing  every  IBS  iteration  in  channe1  -ss 

SPOW  =  (1.-ASTEP)aSPOW+ASTEP*^v’:ABS(R(NPOW)))**2)/2 

The  measured  signal  power  is  recursively  computed  as 

SPOW  =  ((NUM-l)*SP0W+CABS(R(NP0W))**2)/NUM 

which  is  mathematically  equivalent  to 


,  NUM 

SP0W  =  NUM  ^  |R!! 

where  the  RT  are  independent  samples  of  the  signal. 


268.  SQ2 


local  variable  pages  (3-14) 

set  initialize  parameters 
used  channel 
used  noise  filter 

SQ2  =  V2. 


269.  SR 

local  variable 
set  &  used  detector 

The  detector  sets  SR  to  be 

SR  =  SIGN  (1.0,  REAL  (ZC)) 


270.  SUM 

local  variable 
set  &  used  ERFC 
set  &  used  Sync 
set  &  used  Semul 

Sync  sets  SUM  as 

SUM  =  SUM  +  CABS(W(I ))**2 

ERFC  sets  SUM  to  be 

SUM  =  0 

SUM  =  SUM  +  A(I)*(T**I)*DEXP(-(X**2)) 
Semul  sets  SUM  as 
SUM  =  0. 

SUM  =  SUM+((AL*ARG)**J)/JFACT 


271.  SUM1 

local  variable 
set  &  used  Noise 

SUM1  is  set  by  Noise  to 

SUM1  =  1 
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272.  SVAR 


\ 

local  variable  pages  (5-4  back,  3-105) 

used  Gauran 

set  &  used  Noise 

Gauss  in  Random  Number  Generator  Variance 


SVAR  is  the  real  variance  computed  by  NOISE  and  used  in  the  call  to 
Gaussian  random  number  generator  GAURAN 

SVAR  =  TSR/(2. BSR*SNR) 


273.  T 

local  variable 
set  &  used  ERFC 

T  =  (1+. 3275911*X) 

T  =  ONE/ (0NE+P*X) 


274.  TEMPF 

local  variable 
set  &  used  detector 

The  detector  sets  TEMPF  as 

TEMPF  =  AIMAG(F) 


92 


275.  TI 


common  block 
set  DFES 

output  update  input 
set  &  output  Sync 
set  Parin 

Timing  Interval 

When  bit  synchronization  is  required  (NA  >  0  and  SMODE^O),  the 
initial  value  of  TI  is  0  and  the  final  TI  value  is  computed  by  the 
subroutine  SYNC.  If  bit  synchronization  is  not  used,  TI  is  an  input 
parameter  to  the  subroutine  PARIN.  TI  is  a  real  variable  with  the  range 

-1  <  TI  <  1 

Once  computed  by  SYNC  or  input  by  PARIN,  TI  remains  constant. 

Sync  sets  TI  as 

TI  =  T0P/(RMAX*2. )  and  TI  =  0 

TI  is  used  to  compute  the  SINK  array  for  the  Interpolator 
subprogram. 

DFES  sets  TI  =  0.0 


pages  (3-8  back,  3-11,  3-17,  3-19, 
3-33,  3-55,  3-58,  3-83) 


276.  TOP 

local  variable 
set  &  used  Sync 

Sync  sets  TOP  as 

TOP  =  RS(IMAX+1)  if  IMAX  =  1 

TOP  =  RS(IMAX-l)  if  IMAX  =  NPN 

TOP  =  RS(IMAX+1)-RS(IMAX-1)  if  IMAX  >1  and  IMAX  <NPN 
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277.  TSR 


common  block  pages  (3-4  back,  3-12,  3-18,  3-23, 

used  initialize  parameters  3-25,  3-29,  3-49,  3-55,  3-57, 

used  differential  decoder  3-58,  3-66) 

set,  used,  &  output  Parin 
used  Noise 

Tape  Sample  Rate  in  Hz 

TSR  is  a  real  input  constant  to  the  PARIN  subroutine.  It  has  a 
default  value  of  192K  and  must  be  an  integer  multiple  of  BSR. 

TSR  is  used  by  the  Initialize  Parameters  subroutine  in  computing 
NTB,  SIGMA,  the  number  of  tape  samples  per  receiver  sample  NTR  =  TSR/RSR, 
and  the  channel  tap  integer  array  KSET(I)=TSR*DELAY(I)+0. 5.  The  fiOISE 
subroutine  uses  it  to  compute  the  variance 

SVAR  =  TSR/(2*BSR*SNR) 

and  _ 

SIGMA  =  VTSR/(2*BSR*SNR) 


278.  V 

array  pages  (3-25  back,  3-12,  3-17,  3-37, 

set  DFES  3-39,  3-40,  3-54,  3-59) 

set,  used,  &  output  noise  filter 
set  &  output  Parin 

Noise  Filter  Array 

V(I),  1=1,2, — MTAP  is  a  complex  array  used  by  the  Noise  Filter 
subprogram  in  computing  the  noise  filter  output  array  X  when  MTAP>1.  It 
has  a  maximum  size  V(25).  V  is  initialized  by  PARIN  where  it  is  re¬ 
defined  as  a  real  array  (V(50).  When  MTAP>1,  V  may  be  optionally 
initialized  as  an  input  parameter  array  to  PARIN.  The  Default  values 
are  V(I)=0,  1=1,2,-— MTAP. 

DFES  sets  V  as 

V(I)  =  (0.0,  0.0) 

The  Noise  filter  sets  V  as 

V(I)  =  V(I )-VSTEP*F*CONJG(CJ) 
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279.  VIN 


virtual  array 

set  6 1  used  Taper  (Read) 

Complex  data  samples  array 

VIN(1, J)  =  FLOAT ( IDATA{ I ) )*SCALE 
VIN(2,  J)  =  FL0AT(IDATA(I+1))*SCALE 


280.  VKAL 

array  pages  (3-25  back,  3-89,  3-91,  3-92, 

set,  used,  &  output  Fwate  3-96,  3-98) 

used  Bfilt 

Kalman  Update  Vector 

VKAL(I),  1=1,2, - NLTAP  is  a  complex  array  computed  by  FWATE  for 

use  in  updating  the  weight  vector  W  and  array  PVAR  when  NALG=2.  BFILT 
also  uses  VKAL  to  update  BETA  when  NALG=2.  VKAL  is  a  calling  sequence 
argument  of  FWATE  and  BFILT  and  has  a  maximum  size  VKAL(50). 

Fwate  sets  VKAL  to 

VKAL(I)  =  (0.0.) 

if  J>NTAP  then 

VKAL(I)  =  VKAL( I )+PVAR( I , J)*GBACK( JSET( J-NTAP) ) 
if  J<NTAP  then 

VKAL(I)  =  VKAL(I)+PVAR(I ,J)*SBAR(J) 


281.  VMAX 

local  variable 
set  DFES 

set  &  used  noise  filter 
DFES  sets  VMAX  to 
VMAX  =0.0 

The  noise  filter  sets  VMAX  to 
VMAX  =  0 

VMAX  =  CMAG  if  CMAG  >  VMAX 
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282.  VR 


array 

used  channel 
set  Taper  (Read) 

Receiver  Sample  Input  Array 

VR  is  the  receiver  sample  array  returned  by  TAPER  each  bit 
iteration.  VR  is  defined  as  a  real  array  dimensioned  CV(2, 100)  where 
( VR1 , I )  =  real  sample  part  and  VR(2,I)  =  imaginary  sample  part.  It  is 
used  to  set  the  transmit  signal  array  for  the  Channel  subprogram 

Q(K)  =  CMPLX(VR(1,K),VR(2,K), 

k  =  1,2,---NTB+NRQ 

When  IBS=1 ,  TAPER  (READ)  returns  NR  samples  in  VR.  For  IBS  >  1  NTB 
samples  are  returned. 

VR(1,NV)  =  0 

VR(2,NV)  =  0 

VR( 1,1)  =  VIN(l.NVIN) 

VR(2 , I )  =  ViN(2 ,NVIN) 

VR( 1 ,NV)  =  VIN(l.NVIN) 

VR(2 ,NV)  =  VIN(2,NVIN) 


283.  VSTEP 

local  variable 

set,  used,  &  output  noise  filter 
Adaptation  algorithm  step  size 

VSTEP  =  2. *PI*NFLB*RSR/(BSR*BRF*NRB) 
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284.  W 

array  pages  (3-26  back,  3-17,  3-19,  3-41, 

set  OFES  3-42,  3-52,  3-53,  3-59,  3-87, 

output  update  input  3-89,  3-92) 

used  forward  filter 

output  differential  decoder 

set  &  used  Sync 

set  &  used  Fwate 

set  Parin 

Forward  Filter  Weight  Vector 

W(I),I=1,2, — NTAP  is  an  input  parameter  array  to  the  subroutine 
PARIN,  with  the  default  W(l)=l,  W(I)=0,  1=2, 3, --NTAP2.  W  has  a  maximum 
dimension  of  25.  The  array  W  is  updated  each  bit  symbol  iteration  by 
FWATE  using  either  the  Kalman  or  LMS  algorithm.  W  is  used  in  computing 
the  forward  filter  output  array  Y. 

Sync  sets  W  to 

W(I)  =  C0NJG( RSYNC(2*IMAX- ISET ( I )/RSYNC( IMAX) ) 

FWATE  sets  W  to 

W(I)  =  CMPLX(RW.QW) 


285.  WMAX 

local  variable 

set  &  used  forward  filter 

The  forward  filter  sets  WMAX  to 

WMAX  =  0. 

WMAX  =  CABS(W(I))  if  CABS(W(I))>WMAX 
WMAX  =  WMAX*RTHOLD 


286.  WMULT 

local  variable 
set  &  used  Fwate 

WMULT  is  set  by  Fwate  to 

WMULT  =  1. 

WMULT  =  1. -DELTA  if  ALGOR  =  4 
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287.  X 


a)  virtual  array  pages  (3-16,  3-37  to  3-40) 

set  &  used  interpolator 
set,  used,  4  output  noise  filter 
used  forward  filter 
set  4  used  ERFC 
used  Sync 
used  Fwate 

Noise  Filter  Output 

X  is  a  complex  array  of  size  NX=NS  with  a  maximum  dimension 
X(500).  For  MTAP=1  it  is  formed  directly  from  the  Interpolator 
array  S.  When  MTAP  >  1,  X  is  computed  by  the  Noise  filter  sub¬ 
program  each  bit  symbol  iteration.  It  is  used  in  forming  the 
forward  filter  output  array  Y. 

X  is  a  calling  sequence  argument  of  the  subroutines  SYNC  and 
FWATE.  In  SYNC  it  is  used  in  forming  the  synchronization  array 
RSYNC  and  FWATE  uses  it  to  compute  the  SBAR  array. 

The  interpolator  sets  X  to 

X(K)  =  X(K-NRB) 

ERFC  sets  X  =  Z 

The  noise  filter  sets  X  to 

X(I)  =  S(I)  where  I  =  1,NX 

X(K)  =  (0. ,0.) 

X(K)  =  X(K)  +  V(I)  *  XFREQ 

X(K)  =  SLK  +  IBDEL)  -  X(K) 


b)  local  variable  pages  (4-6  back,  3-70,  3-109) 

set  4  used  Gauran 
set  Opgen 
used  Sine 

HSQ  SINC  Input 

Gauran  sets  X  to 

X  =  SQRT(-2.0*SVAR*ALOG(A)) 

Opgen  sets  X  to 

X  =  CMPLX(D(1) ,D(2)) 

X  is  the  real  variable  input  to  the  function  SINC 
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local  variable 
used  channel 


pages  (4-7  back,  3-27) 


XC  corresponds  to  S  in  the  call  to  Gauran 
Channel  Gaussian  Random  Number 

XC  is  the  complex  Gaussian  random  number  returned  by  GAURAN  to  the 
Channel  subprogram.  It  is  used  in  updating  the  complex  H  array. 


289.  XERR 

local  variable 

set  &  used  differential  decoder 

The  differential  decoder  sets  XERR  to 

XERR=AERR*XERR+(1. -AERR)*DR  if  REF=3 

290.  XF1 

local  variable 

set  &  used  noise  filter 

291.  XF2 

local  variable 

set  &  used  noise  filter 


292.  XFREQ 

local  variable 

set,  used,  &  output  noise  filter 
XFREQ  =  (0. ,0. ) 

XFREQ  =  AIX*XF1-CIX*XF2+BIX*S(K) 


293.  XKSET 

local  variable 
set  &  used  Parin 

XKSET  is  set  by  Parin  to 


XKSET  =  TSR*DELAY(I)+0.5 


294.  XX 


local  variable  pages  (4-7  back) 

set  &  used  update  input  _ 

Interpolator  SINC  Function  Input 

XX  is  the  real  variable  input  to  the  function  SINC  for 
computing  the  Interpolator  SINK  array 

XX  =  FLOAT(I)/NTR-TI 


295.  Y 

array  pages  (3-41,  3-43,  3-44) 

set  &  used  forward  filter 
used  &  output  compressor 

Y(K) ,  Input  array  (Complex,  K=1,2,...NY,  NY=NRB  <  80. 
The  forward  filter  sets  Y  to 
Y(K)  =  (0. ,  0.) 

Y(K)  =  Y(K)  +  W(I)  *  S(J  +  IBDEL) 

Y(K)  =  Y(K)  +  W(I)  *  X(J) 

Y(K)  =  Y(K)  -  X(K) 


296.  YC 

local  variable 
set  &  used  channel 

YC  =  HP(J) 
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297.  2 


local  variable  pages  (3-44  to  3-46,  3-50) 

set  &  used  compressor 

output  differential  decoder 

used  ERFC 

used  detector 

Z  =  COMPRESSOR  output  (Complex). 

The  cross  correlation  operation  which  accomplishes  the  bandwidth 
compression  is 

-1  NY 

2  =  (NY)  1  1  PN( I )*Y( I ) . 

1=1 

Z  is  initially  set 

Z  =  (0. ,0)  by  the  compressor. 


298.  ZC 

local  variable  pages  (3-52) 

set  &  used  detector 

output  differential  decoder 

Predecision  Sample  of  Detector 

The  predecision  sample  is  the  sum  of  the  compressor  and  backward 
filter  outputs,  viz.,  as  set  in  the  Detector  subprogram 

ZC  =  KGAIN*Z+C 


299.  ZERO 

algorithm 

local  variable 
set  &  used  Parin 

Parin  sets  ZERO  using  a  Data  Statement  to 
ZERO  =0.0 
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CROSS-REFERENCE  LIST 


KEY 

A  =  set,  used,  and  output 
B  =  set  and  used 
C  =  used  and  output 
D  =  set  and  output 
0  =  output  (printed  out) 

S  =  set 
U  =  used 
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